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1  INTRODUCTION 

As  part  of  a  continuing  effort  to  supply  the  user  community  with 
the  best  phenomenology  codes  available  for  the  evaluation  of  nuclear  ef¬ 
fects  on  sensors,  DNA  has  undertaken  the  task  of  expanding  the  ROSCOE 
(Radar  and  Optical  Systems  Code  with  Nuclear  Effects)  computer  program"*" 
to  include  an  optical/SWIR  capability. 

An  initial  program  planning  phase  was  initiated  in  1976  to  formu¬ 
late  the  objectives  of  this  model,  place  bounds  on  the  scope  and  level 
of  detail  required,  and  suggest  a  plan  for  its  development.  The  plan 
called  for  code  development  during  1977  and  1978,  with  a  completion  date 
of  approximately  January  1979,  followed  by  a  code  validation  and  documen¬ 
tation  effort. 

General  Research  Corporation,  as  part  of  a  team  of  contractors,  has 
participated  in  this  program  development  effort.  GRC's  responsibilities 
have  been  primarily  concerned  with  the  development  of  the  optics  system 
model  and  the  interfacing  of  the  optics  package  (both  systems  and  physics) 
with  the  code  as  a  whole. 

This  paper  presents  a  summary  of  the  optics  code  structure.  In 
Section  2,  a  description  of  the  optical  code  interface  structure  is  given; 
Section  3  presents  a  review  of  one  of  the  major  tasks  performed  by  GRC — 
the  modeling  of  inhomogeneous  regions.  Appendix  A  provides  some  analyti¬ 
cal  backup  material  for  Section  3,  and  Appendix  B  contains  flow  diagrams 
for  the  major  programs  and  subroutines  which  comprise  the  optics  code 
interface  structure. 
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2  SUMMARY  OF  CODE  DEVELOPMENT 

The  following  sections  describe  the  optics  code  structure,  the 
system  simulation  methodology,  and  some  sample  outputs. 

2.1  OPTICS  CODE  STRUCTURE 

The  ROSCOE  code  is  divided  into  a  number  of  separate  computational 
blocks  which  perform  various  parts  of  the  systems  or  physics  simulation. 
Each  block  is  made  up  of  a  number  of  subroutines,  and  is  stored  in  a 
separate  overlay.  These  blocks  are  referred  to  as  "events."  In  the 
execution  of  the  code,  a  series  of  computations,  or  events,  is  set  up. 
Each  of  these  events  is  given  a  calculation  time,  and  they  are  processed 
(i.e.,  the  proper  overlay  is  called  and  executed)  in  a  time-ordered 
fashion.  An  initial  list  of  events  is  set  up  in  the  input  deck  by  the 
user.  Additional  events  may  be  added  to  the  list  as  computations  are 
carried  out  to  complete  the  simulation  specified. 

A  block  diagram  showing  the  major  computational  blocks  and  infor- 

* 

mation  flow  for  the  optics  code  is  shown  in  Figs.  1  and  2.  The  optics 
systems  blocks,  or  events,  are  shown  in  Fig.  1,  and  the  physics  events 
in  Fig.  2.  Some  interface  blocks  which  feed  phenomenology  data  to  the 
systems  blocks  are  shown  in  both  figures. 

The  physics  portion  of  the  model  (Fig.  2)  is  comprised  of  a  burst 
event  which  creates  the  initial  burst  parameters  and  several  update 
events  to  update  the  nuclear  environment  in  time.  In  the  low-altitude 
region  (<90  km),  an  update  event  is  used  to  maintain  the  fireball  and 
debris  region  data  at  two  times.  The  update  event  calls  the  fireball 
and  debris  model  routines  at  specified  times  to  update  these  data,  and 
then  they  are  interpolated  at  intermediate  times  as  requested  for  systems 
calculations  or  for  printout. 


A  similar  diagram  for  the  ROSCOE  radar  code  is  given  in  Volume  1  of  The 
ROSCOE  Manual. 


The  high-altitude  data  are  also  maintained  at  two  times.  Here, 
because  the  nuclear  effects  may  be  more  widespread,  a  grid  of  physical 
data  (chemistry,  striations)  is  computed  in  addition  to  fireball  and 
debris  properties  (which  are  treated  as  overlays  on  the  grid). 

The  optics  system  portion  of  the  code  consists  of  (1)  an  attack 
generation  event,  which  initializes  the  ambient  atmosphere  and  magnetic 
field  models  and  creates  objects  to  be  viewed  by  the  sensor;  (2)  an 
optics  look  event  (OLOOK),  which  sets  up  the  sensor  look  geometry;  (3) 
an  optics  propagation  path  event  (OPROP) ,  which  sets  up  the  number  and 
position  of  sight  paths  in  the  sensor  field  of  view  for  creating  the 
"scene"  at  the  focal  plane;  (4)  a  sight  path  calculation  event  (SPCALC), 
which  performs  the  integration  of  the  radiation  transport  equation  along 
each  sight  path;  (5)  a  simulated  optics  event  (OPSIM),  which  generates  a 
focal  plane  array  using  the  sight  path  data  and  simulates  the  scanning 
and  signal  processing  functions  of  the  optical  sensor;  and  (6)  a  set  of 
system  events  (currently  there  is  only  an  optical  track  function  (OTRACK) 
modeled  here)  to  simulate  various  system  functions  such  as  track,  homing, 
or  discrimination. 

The  systems/physics  interface  structure  is  much  the  same  as  for  the 
radar  code.  The  propagation  and  sight  path  calculation  events  access 
phenomenology  region  data  (e.g.,  fireball,  debris,  dust,  beta  tube  para¬ 
meters)  by  calling  a  subroutine  which  interpolates  these  data  in  time. 

The  sight  path  calculation  event  also  accesses  a  point  properties  module. 
It  supplies  the  time  and  location  of  the  point  in  space,  and  the  point 
properties  module  returns  the  appropriate  data  requested.  These  data 
include  temperature,  mass  density,  chemical  composition,  cloud  statistics, 
striation  statistics,  earth  albedo,  prompt  and  delayed  ionization,  and 
particulate  optical  properties.  In  addition,  the  sight  path  calculator 
accesses  preprocessed  band  emission  data  from  external  data  blocks. 


2.2  SYSTEMS  SIMULATION  METHODOLOGY  AND  SAMPLE  OUTPUT 


This  section  describes  the  optical  system  simulation  methodology. 

The  discussion  which  follows  describes  the  use  of  the  SWIR  code  by  guid¬ 
ing  the  reader  through  the  simulation  procedure  block  by  block,  and 
displaying  some  intermediate  code  outputs. 

1.  Surveillance  Sensor  Dynamics  (OLOOK) .  An  example  of  how 

the  optics  look  (OLOOK)  event  sets  up  the  look  geometry  for  a  surveillance- 
type  sensor  is  illustrated  in  Fig.  3.  The  boresight  of  the  sensor  is 
determined  from  the  sacellite  position  and  a  reference  position,  as  shown. 
The  frame  time  (the  time  it  takes  for  the  sensor  to  complete  its  circular 
scan)  is  input,  and  must  be  consistent  with  the  scan  rate  desired. 

Scanning  motion  is  simulated  by  creating  a  frozen  frame  of  the 
environment  at  a  selected  engagement  time  and  then  sweeping  over  this 
frame  in  a  circular  fashion  at  the  appropriate  rate. 

2.  Object  Dynamics  (QPROP).  In  the  optics  propagation  (OPROP) 
event,  the  object  motion  in  the  sensor  field  of  view  is  treated  by  inter¬ 
polating  in  time  for  the  position  of  each  object,  and  then  transforming 
these  coordinates  to  the  sensor  coordinate  frame.  Objects  treated  include 
nuclear  bursts  and  their  associated  ionization  regions,  booster  plumes,  and 
natural  clouds  (see  Fig.  4) . 


The  sensor  coordinate  frame'  is  defined  by  setting  the  z-vector 
anti-parallel  to  the  boresight.  The  y-vector  is  in  the  plane  formed  by 
the  z-vector  and  the  satellite  position  vector,  and  is  directed  as  close 
to  the  satellite  position  vector  as  possible.  The  x-vector  is  chosen  to 
complete  the  right-hand  system. 


The  block  names  are  given  here  so  that  the  radar  can  refer  back  to  the 
block  diagrams  shown  in  Figs.  1  and  2. 

This  coordinate  system  is  used  in  generating  the  focal  plane  array. 
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•  BORESIGHT  ESTABLISHED  BY  POINTING  AT  REFERENCE  POSITION 

•  FRAME  TIME  BASED  ON  ROTATION  RATE 

•  CIRCULAR  SCAN  SIMULATES  SENSOR  MOTION 
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Figure  3.  Surveillance  Sensor  Dynamics  (OLOOK) 
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•  OBJECTS  INCLUDE  -  NUCLEAR  BURSTS  (ASSOCIATED  REGIONS) 

-  BOOSTER  PLUMES 

-  NATURAL  CLOUDS 

•  SENSOR  COORDINATE  FRAME 


Figure  4.  Object  Dynamics  (OPROP) 
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3.  Object  Characterization  (OPROP) .  To  characterize  the  objects 
in  the  focal  plane  of  the  sensor,  the  code  integrates  the  radiation  trans¬ 
port  along  a  number  of  paths  from  the  sensor  which  intersect  each  object. 
Since  these  path  integrations  can  be  time-consuming,  an  attempt  is  made  to 
minimize  the  number  of  paths  used  by  selecting  a  relatively  few  paths  and 
interpolating  between  them  later  to  generate  finer  detail  in  the  focal 
plane. 

Objects  or  regions  are  divided  into  three  types:  (1)  point  objects, 
which  are  objects  whose  dimensions  are  small  compared  to  the  sensor  reso¬ 
lution;  (2)  extended  objects  (such  as  fireballs),  which  may  be  bright  and 
cover  many  resolution  cells;  and  (3)  background,  or  other  regions  which 
are  less  bright  or  very  extensive  (cover  the  entire  field  of  view). 

Point  objects  can  be  characterized  by  a  single  path,  as  shown  in 
Fig.  5.  Extended  objects  are  treated  by  taking  the  object  projection  in 
a  plane  normal  to  the  vector  between  the  sensor  and  the  object,  and  select¬ 
ing  a  set  of  points  within  this  boundary.  Finally,  the  background  is 
treated  in  as  much  (or  as  little)  detail  as  the  user  feels  is  warranted 
by  setting  up  a  grid  of  points  within  the  sensor  field  of  view.  The  user 
selects  the  grid  spacing  in  the  input  deck.  A  default  value  of  four  points 
is  used  as  a  minimum. 

4.  Number  of  Paths/Object  (OPROP).  The  number  of  paths  used  to 
characterize  an  extended  object  is  a  function  of  the  object  size  and  the 
sensor  resolution,  as  shown  in  Fig.  6.  One,  five,  seven,  or  thirteen 
paths  through  the  region  are  set  up  for  integration  using  the  scaling 
rules  shown. 


5.  Example  Sight  Paths  to  a  Fireball  (OPROP) .  An  example  set  of 

sight  paths  used  in  generating  the  focal  plane  image  of  a  fireball  is 
shown  in  Fig.  7.  Here  a  striated  fireball  is  viewed  by  a  sensor  on  a 
synchronous  satellite.  The  sensor  resolution  was  small  compared  to  the 
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NPX  =  IFIX  (aX/ RESOL)  -  1 
NPY  =  IFIX  (aY/RESOL)  -  I 
WHERE:  AX, AY  =  OBJECT  ANGULAR  SUBTENSE 
RESOL  =  SENSOR  ANGULAR  RESOLUTION 

y 


1  POINT 


5  POINTS 


7  POINTS 


NPX  =  1 
NPY  =  1 


2  <  NPX  <  4 
2  <  NPY  <  4 


2  <  NPY  <  4 


13  POINTS 


4  <  NPX 
•X  4  <  NPY 


Figure  6.  Number  of  Paths /Object 
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Figure  7.  Example  Sight  Paths  for  Fireball 
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object  dimensions,  so  thirteen  paths  were  used  (denoted  by  dots  in  the 
figure) .  A  symmetric  set  of  paths  were  set  up  in  this  case  to  model  the 
spherical  region. 


6.  Path  Integration  (SPCALC) .  As  mentioned,  for  each  sight  path 

specified  within  the  field  of  view,  the  equation  of  radiation  transport 

must  be  integrated  from  the  satellite  to  the  ground.  The  methodology  de- 
2 

ve.loped  in  WOE  II  for  optical  sight  path  integration  through  multiple 
overlapping  regions  of  disturbance  has  been  adopted  for  use  in  ROSCOE. 


The  methodology  involves  generating  a  cell  list  of  integration  1 

regions  along  the  path  (Fig.  8) .  These  regions  may  be  further  subdivided  f 

into  integration  elements.  The  cell  list  consists  of  bounding  cells  at 
the  sensor  and  at  some  terminator  (the  ground,  in  this  case),  atmospheric  | 

shells  of  about  a  scale  height  in  width,  and  disturbed  region  cells.  A  f 

ranking  scheme  for  determining  the  relative  importance  of  overlapping  | 

regions  has  also  been  taken  from  WOE.  ! 

j 

7.  Sight  Path  List  Cell  Types  (SPCALC).  Table  1  shows  the  • 

various  sight  path  cell  types  that  may  be  generated  in  the  code  and  the  t 

logical  rules  used  in  overlapping  regions.  Each  cell  type  has  a  region  ‘ 

I 

code  and  is  given  a  rank;  the  higher  the  rank,  the  more  important  the 
region.  The  region  codes  for  all  nuclear  disturbed  regions  are  identified 
by  the  number  1000  and  distinguished  by  adding  a  region  index.  The  region 
index  is  merely  a  counter  of  the  number  of  "region  kinds"  intersected. 

For  example,  the  first  disturbed  region  intersection  has  an  index  of  1  and 
a  region  code  of  1001,  and  so  forth.  : 


8.  Example  Cell  List  (SPCALC).  An  example  cell  list  generated 
in  the  code  for  a  typi.cal  sight  path  is  shown  in  Fig.  9.  The  following 
information  is  given: 
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• . . . . . . * . . 


•  CELL  LISTS  CONSIST  OF  —  (1)  BOUNDING  CELLS 

(2)  ATMOSPHERIC  SHELLS 

(3)  DISTURBED  REGION  CELLS 

•  RANKING  SCHEME  USED  FOR  OVERLAPPING  REGIONS 


Figure  8.  Path  Integration  (SPCALC) 


TABLE  1 

SIGHT-PATH  LIST  CELL  TYPES 


Sljht  Path  Region  Region  Coce  t*t^RAHK* 


Ceosetric  Overlap 
Logical  Suits  lK£tPs 


Cloud  layer 


None  of  lower  RANK 

Hone  of  lower  RANK 

None  of  lower  RANK 

None  of  tower  RANK 

None  of  lower  RANK 

All  regions  of  lower  RANK 
except  aobient  atmosphere 


Ambient  atmosphere 

Ambient  chemiluminescent 
pancake 


Ion-leak  and  CKiX  pitch 

Low-altitude  cheailusinesceit 
spr.tre 

High-altitude  cheailuaines- 
cent  pancake 

Beta  tube 


Beta  sheath 


Shock  heated  vortex 


Dust  cloud  stem 


Dust  cloud  cap  and  condensed 
debris  particulates 

Spheroidal  weapon  debris, 
vapor  and  particulates 

Toroidal  weapon  debris, 
vapor  and  particulates 

Fireball  wake 

Ocbris  loss-cone  patch 


Shell  index 

NDS  ♦  1001 

Region  index 
♦  1000 


11  to  IS 


All  regions  of  lower  RANK 
Debris  particulates 


None  of  lower  RANK 


Debris  particulates 
Debris  particulates 


Debris  particulates 

Debris  particulates,  fire¬ 
ball  wake,  and  all  other 
beta  regions 

Ocbris  particulates,  fire¬ 
ball  wake,  and  all  other 
beta  regions 

Debris  particulates,  chemi¬ 
luminescent  regions,  fire¬ 
ball  wake,  and  beta  regions 

All  other  dust  and  particu¬ 
late  regions 

All  other  dust  and  particu¬ 
late  regions 

All  regions  of  lower  RANK 
except  aebient  atcasphere 

All  regions  of  lower  RANK 
except  ambient  atcasphere 

All  regions  of  lower  RANK 

Debris  particulates 


*lho  RANK  of  disturbed  regions  decreases  linearly  with  age  to  a  value  that  is  almost  one  below  the 
Initial  value  when  the  age  reaches  One  hour  ane  thereafter  remains  constant. 

^Append  to  all  entries  ’and  regions  of  greater  RAN K.' 

clhe  chemiluminescent  pancalc  is  giver,  a  fixed  age  of  ten  minutes. 
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Figure  9.  Example  Cell  List 


'■  *'“M.  >>  v;n>i  i  /!.•■.,  iytr  ii'fHWSP’WW'Wft'M  t.’rw  ijjjiv ,  Vn,t. 


I  Cell  number 

LINK  A  packed  word  used  to  find  addresses  of  the  previous 

and  next  cells  (datasets) 

DSPWRD  Pointer  word  in  the  DSA  scheme  which  allows  one  to 

access  the  dataset  for  this  cell 

ICREG  Region  code  for  sight  path  cell  terminating  on  this 

boundary 

DISTIC  Distance  alcng  the  sight  path  from  the  sensor  to  this 

cell  boundary 

KEEPR  A  packed  word  containing  subordinate  region  indices 

(up  to  five  subordinate  region  indices  can  be  stored) 

KIND  (last  column  shown)  Disturbed  region  KIND 

The  cell  list  contains  null  cells  (ICREG  =  0)  at  the  sensor  and  at 
the  ground  (DISTIC  =  3.5824E9)  and  twenty  atmospheric  cells  (ICREG  =  1-20) 
from  the  ground  to  100  km  (DISTIC  =  3.5724E9).  A  target  cell  (ICREG  = 
-50)  of  zero  length  has  been  inserted  at  100  km,  and  six  disturbed  region 
cells  (ICREG  >  1000)  have  been  added  between  66  km  and  75  km. 


The  cell  list  would  then  be  used  as  a  guide  for  the  integration  of 
the  radiation  transport  along  this  path.  Each  path  would  be  integrated 
and  the  resulting  information  passed  to  the  focal  plane  generation  block 
of  the  code. 

9.  Focal  Plane  Generation  (OPSIM).  To  create  a  focal  plane 
image  of  the  environment,  the  code  grids  the  field  of  view  into  a  fine 
mesh  (the  mesh  size  is  taken  as  something  smaller  than  a  detector  diameter), 
and  then  uses  the  sight  path  data  computed  at  selected  points  to  map  the 
intensity.  Recall  that  we  previously  mentioned  that  objects  or  regions 
in  the  field  of  view  are  divided  into  three  types:  background  (grid), 
extended  objects,  or  point  objects,  and  that  each  type  was  characterized 


by  one  or  more  sight  paths.  These  are  used  as  follows  in  the  focal  plane 
development  (Fig.  10): 

a.  The  background  is  generated  by  taking  the  sight  paths  asso¬ 
ciated  with  the  grid  and  using  bivariate  linear  interpolation 
to  compute  the  values  at  each  cell  in  the  mesh. 

b.  Ext 'Med  objects  are  then  added  one  at  a  time  by  taking  the 
sight  paths  associated  with  each  object  and  following  this 
procedure: 

•  A  rotated,  magnetic,  field-aligned  coordinate  system  is 
set  up,  called  the  F-P  frame - 

•  The  maximum  and  minimum  region  dimensions  in  the  F-P 
frame  are  determined. 

•  For  each  cell  in  the  F-P  grid,  the  code  determines  if 
the  sight  path  through  this  point  intersects  the  object. 
If  it  does,  the  three  path  datasets  nearest  this  point 
are  found  and  used  to  interpolate  for  the  mean,  the 
standard  deviation,  and  the  correlation  lengths  of  the 
intensity  distribution  at  the  point.  A  normalized  random 
sample  is  generated  using  the  correlation  lengths  and  a 
one-point  recursion  formula.  The  intensity  is  then  com¬ 
puted  as  the  mean  plus  the  standard  deviation  times  the 
random  sample. 

•  The  point  is  transformed  back  to  the  sensor  coordinate 
frame  and  the  intensity  at  this  point  is  taken  as  the 
maximum  of  the  computed  value  and  the  old  array  value. 

c.  Finally,  point  objects  are  added  to  the  array  at  the  appro¬ 
priate  locations  specified  by  their  sight  path  location  in 
the  field  of  view. 
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BACKGROUND 


•  PLANAR  INTERPOLATION  OF  GRID  DATA 


OBJECTS  —  MEAN 


•  OBJECTS  ADDED  ONE  AT  A  TIME 

•  PLANAR  INTERPOLATION  OF  OBJECT 
DATA  FOR  MEAN  INTENSITY 

•  CONSTRUCTED  IN  COORDINATE  SYSTEM 
ALIGNED  WITH  FIELD 


—  STRUCTURE 


POINTS 


•  PLANAR  INTERPOLATION  FOR  SIGMA 
AND  CORRELATIONS 

■  RANDOM  SAMPLES  GENERATED  USING 

ONE-POINT  RECURSION  AND  CORRELATIONS 

0  INTENSITY  AT  POINT  FOUND  BY  COMBINING 
MEAN  AND  RANDOM  COMPONENTS 

0  TRANSFORMATION  BACK  TO  SENSOR 
FRAME 

0  FINAL  VALUE  AT  GRID  POINT  TAKEN 
AS  MAX  OF  CURRENT  AND  OLD  VALUES 


ADDED  AT  APPROPRIATE  LOCATIONS 
IN  ARRAY 


Figure  10.  Focal  Plane  Generation  (OPSIM) 
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10.  Examples — Focal  Plane  Array  (OPSIM).  A  number  of  test  runs 

were  conducted  to  determine  whether  the  code  structure  up  to  this  point 

produced  reasonable  outputs.  The  procedure  used  in  this  verification 

process  consisted  of  (1)  inputting  fictitious  values  for  the  mean  and 

* 

variation  of  the  integrated  intensity  for  each  sight  path,  (2)  using 
these  data  to  generate  the  focal  plane  array,  and  (3)  producing  contour 
plots  of  the  array.  The  contour  plots  are  produced  by  plotting  an  alpha¬ 
numeric  character  or  blank  for  successive  increments  of  intensity  ("A" 
being  the  lowest  intensity,  and  "Z"  the  highest).  Thus  in  a  region  of 
increasing  intensity,  one  would  see  a  character  (or  a  repetition  of  the 
same  character) ,  followed  by  a  blank  region,  followed  by  a  repetition  of 
a  higher-order  character,  and  so  forth.  A  number  of  these  example  out¬ 
puts  follow. 


Variation  in  Mean  Intensity  (Fig.  11) .  Contour  maps  were  produced 
for  a  case  where  a  single  fireball  was  being  viewed  by  the  sensor.  Two 
sets  of  inputs  were  assumed  for  the  fireball  intensity;  the  first  assumed 
the  mean  was  constant  and  the  variation  was  about  10  percent  of  the  mean, 
and  the  second  assumed  that  the  mean  varied  with  the  depth  of  the  region 
intersected,  and  again  the  variation  was  taken  as  10  percent  of  the  mean. 
Several  key  points  are  shown  in  these  plots: 

*  A  reasonable  representation  of  the  fireball  projection  has 
been  produced.  In  this  case,  the  fireball  dimensions  are 
large  with  respect  to  the  sensor  resolution,  so  the  maximum 
number  of  paths  (13)  were  used. 

♦  The  striated  nature  of  the  fireball  is  easily  seen  in  the 
first  plot.  The  field  direction  is  about  -103°  ccw  from 
the  horizontal  x-axis,  and  correlation  distances  of  0.6  km 
(equivalent  to  0.0167  mrad  in  this  example)  across  the 
field  and  1000  km  along  the  field  were  assumed.  This  shows 
that  the  one-point  recursion  formula  produces  reasonable 
outputs  in  the  two  spatial  dimensions. 
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B)  MEAN  VARIES  WITH  SQUARE  OF  THE  DEPTH  OF  THE  REGION  INTERSECTED 


Figure  11.  Variation  in  Hean  Intensity 
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•  The  striations  are  somewhat  less  distinctive  in  the  second 
case  where  the  mean  falls  off  from  the  fireball  center;  but 
the  representation  appears  to  be  reasonable. 

Spatial  Correlation  (Fig.  12).  The  correlation  distance  across  the 
field  was  varied  to  further  test  the  recursion  formula.  The  plots  show 
the  same  fireball  with  assumed  correlation  distances  of  0.06  km  and  0.6 
km,  respectively.  As  one  could  expect,  the  plots  show  that  as  the  corre¬ 
lation  distance  increases,  the  spacing  between  striations  in  the  realiza¬ 
tion  increases. 

Time  Correlation  (Figs .  13  and  14) .  To  test  the  third  dimension  of 
the  recursion  formula  (i.e.,  time  correlation),  contour  plots  at  look 
times  of  0  and  1  second  were  made  with  two  assumed  correlation  times. 

The  first  set  (Fig.  13)  shows  contours  for  a  case  where  the  corre¬ 
lation  time  (10  s)  is  large  in  comparison  to  the  look  interval  (1  s) .  As 
one  would  expect,  the  structure  within  the  fireball  cl  anges  only  slightly. 

The  next  set  (Fig.  14)  shows  contours  for  a  correlation  time  much 
smaller  (0.1  s)  than  the  look  interval.  Here,  the  random  samples  generated 
at  the  two  look  times  are  nearly  independent,  so  that  the  realizations  of 
the  structure  are  completely  different. 

Multiple  Fireballs  (Fig.  15).  Next,  a  run  was  performed  to  test 
the  logic  for  generating  overlapping  fireballs.  In  the  upper  part  of  the 
figure,  the  single  fireball  used  previously  is  shown.  In  the  lower  figure, 
a  second  fireball  of  the  same  size  but  displaced  to  the  side  has  been 
generated.  This  second  fireball  has  the  same  mean  intensity  and  structure. 

The  plot  shows  that  the  fireball  profiles  are  clearly  defined,  and 
a  reasonable  realization  of  the  structure  has  been  generated. 
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3)  TWO  OVERLAPPING  FIREBALLS 

Figure  15.  Multiple  Fireballs 
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Multiple  Objects — growth  with  time  (Figs.  16  and  17).  Another  test 
of  the  multiple  object  logic  was  made  where  two  different  object  types 
were  assumed:  a  fireball  and  a  beta  tube.  For  simpJicity,  the  fireball 
and  beta  tube  were  assumed  to  have  the  same  mean  intensity  and  intensity 
variation  due  to  structure.  The  plots  (Fig.  16)  show  (1)  the  fireball 
alone,  and  (2)  the  fireball  and  beta  tube  together  at  a  time  10  seconds 
after  burst. 

A  second  set  of  plots  are  shown  in  Fig.  17  at  a  later  time  (20 
seconds  after  burst)  to  show  that  the  program  can  produce  reasonable 
object  representations  as  the  regions  grow  with  time.  Note  here  that 
both  the  fireball  and  beta  tube  have  expanded  radially  about  30  percent. 

11.  Sensor  Scan  (SPIRE) .  As  mentioned  above,  once  a  realization 
of  the  environment  has  been  created  in  the  focal  plane,  the  sensor  scan¬ 
ning  motion  is  simulated  by  sweeping  over  this  frame  with  a  detector 
array.  The  scanning  motion  and  number  and  position  of  the  detectors  in 
the  array  are  specified  by  the  user.  In  the  example  depicted  in  Fig.  18, 
a  surveillance  sensor  which  uses  a  circular  scan  and  a  simple  linear 
alignment  of  four  detectors  is  assumed.  Because  the  focal  plane  array 

is  gridded  with  a  spacing  on  the  order  of  the  sensor  resolution,  the 
scanning  process  merely  consists  of  picking  up  focal  plane  array  values 
in  the  proper  order  for  each  detector.  The  output  consists  of  a  data 
stream  for  each  detector. 

12.  Optics  Data  Processing  (SPIRE) .  The  sensor  data  processing 
model  consists  of  over  thirty  processing  algorithms  which  can  be  used  in 
a  building-block  approach  for  modeling  a  specific  data  processing  system. 

A  simple  example  of  this  approach  is  shown  in  Fig.  19  for  the  sur¬ 
veillance  sensor  mentioned  above.  For  each  detector,  the  focal  plane 
array  is  scanned  to  produce  a  data  stream  of  samples,  detector  noise  is 
added  co  each  sample,  the  data  stream  is  smoothed,  adjacent  samples  are 
differenced,  and  a  threshold  test  is  applied  to  detect  targets. 
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Figure  17.  Multiple  Objects — Growth  With  Time 
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Figure  18.  Sensor  Scan 
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OUTPUT  A  —  RAW  SIGNAL  AT  SENSOR 
OUTPUT  B  —  NORMALIZED  SIGNAL 


OUTPUT  C  —  SMOOTHED  SIGNAL 


OUTPUT  D  -  DIFFERENCED  SIGNAL 


OUTPUT  E  -  TARGET  DETECTIONS 


Figure  19.  Optics  Data  Processing 


Outputs  are  produced  at  several  stages  in  this  process.  The  raw 
signal  as  it  is  sampled  from  the  focal  plane  is  output,  as  well  as  the 
same  signal  normalized  to  the  sensor  NEFD,  the  signal  after  the  addition 
of  detector  noise  and  smoothing  has  been  performed,  the  differenced  signal, 
and  finally,  a  flag  to  denote  target  detections  (samples  that  exceed  an 
input  threshold) . 

13.  Surveillance  Sensor — Example  Case.  The  sensor  scanning  and 
processing  functions  were  exercised  for  the  following  scenario: 

A  surveillance  sensor  attempts  to  detect  a  booster  which 
has  been  launched  just  prior  to  the  sensor  scan  of  the 
launch  site.  A  nuclear  burst  which  has  occurred  just  10 
seconds  earlier  is  at  an  altitude  just  below  the  booster, 
and  produces  a  high  level  of  background  radiation. 

The  scenario  is  shown  in  Fig.  20,  along  with  the  focal  plane  gen¬ 
erated  by  the  program.  The  fireball  and  beta  tube  were  assumed  to  have 

* 

constant  mean  intensities  of  20  dB  and  intensity  variations  due  to 
striations  of  10  dB.  The  booster  consists  of  a  single  point  in  the  plane 
(the  character  "D"  at  the  center  of  the  plot)  with  an  intensity  of  30  dB. 

The  data  stream  produced  by  the  scanning  and  processing  of  the 
focal  plane  data  is  shown  in  Fig.  21  for  one  detector.  The  output  con¬ 
sists  of  the  time  at  which  each  sample  is  taken,  the  detector  number,  the 
central  wavelength  for  this  measurement,  the  position  of  the  sampled 
point  in  azimuth  and  elevation  of f-boresight,  and  the  five  processing 
outputs  mentioned  above. 

The  data  show  that  the  sensor  is  able  to  detect  the  target  in  the 
presence  of  the  noisy  background  created  by  the  burst. 


The  intensity  contours  are  plotted  in  terms  of  signal/noise.  Signal/ 
noise  is  defined  here  as  the  irradiance/NEFD. 
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Figure  21.  Surveillance  Sensor — Example  Output 
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MODELING  OF  INHOMOGENEOUS  REGIONS 


3.1  INTRODUCTION 

In  the  ROSCOE  simulation  currently  available  for  radar  and  satel¬ 
lite  communications,  two  forms  of  inhomogeneities  are  considered  for  the 
ionized  regions  that  affect  propagation  at  microwave  frequencies.  These 
are  the  fireball  surface  bumpiness  responsible  for  radar  clutter  at  low 

altitudes,  and  field-aligned  striations  within  and  adjacent  to  high- 
A 

altitude  fireballs. 

Extension  of  the  ROSCOE  simulation  to  include  effects  of  nuclear 
detonations  on  infrared  systems  has  necessitated  some  additions  to  the 
set  of  structured  regions  modeled.  Fireballs  at  all  altitudes  are  ob¬ 
served  to  show  flocculent  structure  inside  the  bright  fireball,  at 
frequencies  in  the  visible  range.  This  structure,  usually  attributed 
to  Taylor  instabilities  during  the  initial  disassembly  of  the  warhead 
debris,  should  be  present  in  the  2-5  pm  region  also,  and  at  high  altitude 
evolves  in  time  to  form  striations  aligned  along  the  geomagnetic  field. 

In  addition  to  this  structure  inside  fireballs,  other  irregular 
regions  present  in  the  visible  and  SWIR  include  beta  patch  irregularities 
and  light  reflected  from  clouds.  The  beta  patch  irregularities  are  pre¬ 
sumably  causally  related  to  corresponding  structure  in  the  fission  debris 
at  higher  altitudes  along  the  geomagnetic  field  lines,  and  will  be 
modeled  in  some  as  yet  unspecified  way  within  the  ROSCOE  simulation. 
Reflection  of  solar  and  bomb-thermal  light  from  fireballs  is  being  ac¬ 
counted  for  separately,  so  this  section  deals  only  with  the  modeling  of 
radiating  structures  at  high  altitude,  the  striations  within  and  near 
high-altitude  fireballs. 


Even  the  surface  lumpiness  at  low  altitudes  is  not  modeled  directly, 
but  only  through  its  effects  in  extending  the  reflected  radar  pulse 
length  beyond  the  extent  of  the  specular  reflection  region. 
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ROSCOE  presently  calculates  the  slow  growth  rates  of  striated 
regions,  assumed  field-aligned  at  all  times,  from  the  gradient-drift 
instability  as  described  in  Volumes  6  and  18  of  The  ROSCOE  Manual.  The 
growth  rate  for  this  instability  is 

a  -  (V.  -  V  )  *  Vi  In  n  seconds  ^  (3.1) 

in  e 

where  -  V  is  the  relative  slip  velocity  between  ions  and  neutrals 
and  the  component  of  the  logarithmic  gradient  of  electron  density  ng 
perpendicular  to  the  magnetic  field  is  used.  Thus  striations  should 
form  from  the  gradient-drift  instability,  wherever  a  gradient  in  electron 
density  and  a  slip  velocity  between  neutrals  and  ions  coexist.  Initially 
high-altitude  fireballs  have  strong  electron  density  gradients  near  their 
edges  (although  ROSCOE  treats  high-altitude  fireballs  as  geometric  over¬ 
lays  with  constant  electron  density  across  planes  of  equal  altitude). 

Thus  fireballs  should  striate  on  the  surface  soon  after  magnetic  contain- 

3 

ment  is  attained.  Other  regions  identified  as  possible  sources  for 
striation  growth  by  this  mechanism  are  the  fireball  core  and  the  fireball 
bottom  after  this  region  has  passed  through  its  apogee  and  falls  back 
into  the  atmosphere. 

Thus  high-altitude  fireball  structure  of  relevance  to  radiance 
variations  comprises  initial  flocculence,  presumably  isotropic  at  early 
times,  but  growing  with  time  along  the  magnetic  field;  surface  striations 
around  the  periphery  of  the  fireball;  and  volume  striations  throughout 
the  body  of  the  fireball.  These  may  all  be  modeled  with  different  para¬ 
meters  as  regards  size  spectrum,  degree  of  structure  (variation  of 
radiance),  and  relative  correlation  distances  along  and  normal  to  the 
geomagnetic  field.  However,  the  modeling  must  be  consonant  with  the 
current  representation  of  generalized  high-altitude  structures,  whose 
form  is  discussed  in  the  next  subsection. 


3.2  GENERATION  OF  RADIANT  INTENSITIES  IN  THE  IDEAL  FOCAL  PLANE 

During  the  initial  structuring  of  the  optics  routines  for  ROSCOE, 
a  particular  concern  was  the  question  of  how  statistical  properties  of 
structured  radiating  regions  should  be  represented  in  the  patch  plane, 
or  ideal  focal  plane,  and  how  specific  realizations  of  how  these  statis¬ 
tics  should  be  generated.  The  two  basic  approaches  considered  were: 

1.  Accumulation  of  statistics  along  a  number  of  paths,  each 
corresponding  to  a  particular  point  in  the  patch  plane. 

The  number  of  paths  must  be  sufficient  to  describe  indi¬ 
vidual  region  boundaries  in  the  plane.  A  specific 
realization  is  then  generated  from  a  recursion  formula. 

2.  Calculation  of  statistics  for  a  number  of  lines  of  sight 
within  each  object  separately.  A  specific  realization 
is  then  generated  for  each  object,  from  a  recursion  for¬ 
mula,  and  superimposed  on  the  patch  plane  along  with  the 
specific  realization  for  other  regions. 

For  the  specific  case  of  sensors  at  synchronous  altitudes,  stria- 
tions  in  regions  ac  100-km  altitude  look  about  like  those  at  200-km 
altitude,  because  of  the  great  distance  to  the  satellite  sensors.  It 
would,  however,  be  dangerous  to  capitalize  on  this  fact  by  designing 
the  code  in  such  a  specific  way  that  later  versions  for  other  sensor 
systems  would  have  to  be  restructured  in  some  major  way. 

We  present  below  a  method  for  generating  correlated  intensity  pat¬ 
terns  in  the  patch  plane,  together  with  the  equations  necessary  for 
extrapolating,  interpolating,  and  joining  different  numerical  sequences 
in  several  dimensions,  when  the  variance  of  the  computed  values  is  per¬ 
mitted  to  be  any  arbitrary  function  of  position,  and  when  the  spatial 
autocorrelation  function  has  two  orthogonal  axes  of  symmetry.  The  only  as- 
sumotions  are  that  the  statistics  are  Gaussian  and  that  the  autocorrelation 


39 


function  is  exponential,  with  a  different  (but  constant)  e-folding  length 
along  the  two  space  axes.  This  situation  is  felt  ?  appropriate  to  the 
case  of  an  infrared  sensor  scanning  in  some  arbiri.;  '.rection  along  high- 
altitude  striations,  where  the  magnetic  field  pro.  ..iS  an  automatic  elonga¬ 
tion  of  the  luminosity  parallel  to  the  field.  Time  decorrelation  due  to 
relative  motion  of  the  line  of  sight  through  tl  .  striations  is  also  assumed 
to  be  exponential  in  time. 


Consider  a  two-dimensional  array  of  intensities  ,  which  we 

take  as  having  already  been  "detrended;"  i.e.,  the  smooth  variation  of 
intensity  as  the  line  of  sight  sweeps  through  regions  of  different  total 
luminosity  has  been  subtracted  out  so  that  we  need  consider  only  fluctua¬ 
tions  about  the  mean  value  of  zero.  As  the  line  of  sight  moves,  the 
fluctuations  of  course  change  also,  so  we  ascribe  a  variance  o-j^j  to 
the  signal,  and  furthermore  we  assume  that  there  exist  two  different 
exponential  correlation  distances  in  the  i-  and  j-directions.  The  total 
intensity  (which  is  assumed  detrended)  and  variance  are  computed  from  the 
known  fireball  parameters  along  a  line  of  sight  designated  as  i,j  ,  or 
along  several  lines  of  sight.  Our  problem  then  is  to  compute  a  "reasonable" 
representation  of  the  expected  intensity  fluctuations  at  other  arbitrary 
points,  and  to  ensure  that  all  points  correlate  properly  with  themselves 
and  each  other.  We  would  like  to  use  a  simple  recursion  relation  relating 
the  next  value  at  i+n,i+m  only  to  j  by  an  equation  resembling  that 
for  stationary  Markov  processes.  We  want  2  computed  value  to  have  the 
following  properties: 


a. 


<S..  > 

l+n,  j+m 


2 

°i+n, j+m 


(3.2) 


where  a.  is  the  variance  of  the  signal. 

r+n,j+m  6 


b.  <s-4^  S,  >  =  FnPm  0.  .  0..  .. 

i+n,3+m  i,j  i,j  i+n,j+m 


(3.3) 


A  discussion  of  the  possible  problems  engendered  by  this  choice  is  given 
in  a  later  section,  with  some  suggestions  for  improvements,  if  desired. 
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Hare, 

pn  =  e-nAx/Lx 
pin  =  e-mAy/L2 

where  and  are  the  correlation  distances  along  the  two  axes, 

and  Ax  and  Ay  are  the  linear  separation  of  two  neighboring  points  in 
the  mesh  to  be  computed.  The  third  property. 


c.  <S.  >  =  <S..  ..  >  =  0 

1,3  l+n, j+m 


(3.4) 


is  automatically  satisfied  by  the  detrending. 


Let  us  try  a  recursion  relation  of  the  form 


3+P^+P1.  =  FnPm  fill  +  J'i  _  F2nP2m  R 

o. ,  . ,  a.  .  1  i,  j 

i+n,j+m  i,j 


(3.5) 


where  R.  .  is  a  random  Gaussian  var.'.ate  of  variance  1,  as  are  the  quan- 
i,3 

tities  S.  ./a.  .  .  This  equation  is  seen  to  satisfy  the  three  conditions 
1,3  1,3 

above.  Furthermore,  the  process  may  be  extended  to  the  values  j+m' 

where  n'  >  n  and  m'  >  m  ,  merely  by  successive  application  of  Eq.  3.5. 


S  S 

i+n*  ^j+m*  _  pn'-n  ptn’-m  i+n,,j+m 

i+n',j+ra’  ai+n,j+m 


+  1/1  -  F2(n'~n)  p2(m’-m)  R 


(3.6) 


This  clearly  satisfies  condition  (a),  and  if  we  compute  the  correlation 
with  the  previously  computed  S_^  ,  it  is  seen  to  obey  (b)  too,  since, 
from  Eq.  3.6, 
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l|,[jl|il  I'll1  mill'1  llrliHU  ini*  Util  M-8-1  k  willUtliUuflirfbl  fcl.il 


<s..  ,  ,  s.  > 

i+n  ,_i+m  1,1 

ai+n' , j+m’ 


_n'-n  -m'-m 

- - -  <S..  •  S.  > 

°i+n, j+m  1+n’J*n  ^ 


=  Fn'-n  Pm’-m  FV  a 

i»j 


(3.7) 


Therefore  the  sequence  may  be  automatically  extrapolated  to  larger  values 
of  n,m  in  steps  of  arbitrary  size.  (The  case  where  either  n'  or  m' 
is  less  than  n  or  m  is  a  special  case  of  the  interpolation  equations 
to  be  discussed  next.)  From  symmetry,  it  is  also  seen  that  the  sequence 
may  also  be  extended  in  the  other  direction  from  the  lowest  values  of 
i,j  .  But  now  the  question  arises  as  to  how  to  interpolate;  i.e.,  having 
computed  from  Eq.  3.5,  what  do  we  do  about,  say,  Si+n  j+m-1  * 

The  mere  fact  that  we  have  computed  the  point  at  i+n,j+m  by  drawing 
from  the  random  variable  R  automatically  sets  a  constraint  on  the  va1ue 
of  j+m-1  »  since  it  must  now  correlate  with  si+n  >  as  well 

as  with  S.  .  . 

Let  us  compute  the  cwo  other  corner  points,  i,j+m  and  i+n,j  , 
which,  by  Eq.  3.5,  are  given  by 


(.S/a).  =  Pm(S/o) .  +  \!1  -  PZm  R. 

*>3  ’  3 


(S/a)..  .  =  F (S/a).  .  + 

i+n»3  r>3 


1  -  F2n  R. 

l 


(3.8) 


(3.9) 


where  we  are  henceforth  using  the  subscripts  only  on  the  quantities 


S.  Jo.  .  =  (S/a).  . 

i,j  x,  j  x,j 


These  two  values  have  unity  variance,  and  clearly  correlate  in  the  ap¬ 
propriate  way  with  S-j^j  but  not  with  in  Eq.  3.5.  The 

reason  is  clear.  Whatever  choice  was  made  of  the  random  variate  R,-  • 

X  5  J 

in  Eq.  3.5  sets  a  constraint  on  the  separate  values  R^  and  R_.  in 
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Eqs.  3.8  and  3.9 — which  constraint  may  be  easily  computed.  Thus,  apply¬ 
ing  (3.8)  and  (3.9)  successively. 


(S/c)..  =  Pm(S/cr) .  .  +  Wl  -  ?2m  r! 

l+n, j+m  i+n,3  1  j 

=  PmFn(S/o).  .  +  Pm  Jl  -  F2n  R.  +  Jl  -  P2m  R!  |  (3.10) 

i,J  '  i  1  j) 

This  clearly  satisfies  conditions  (a)  and  (b)  as  applied  to  all  pairs  of 
points.  The  term  in  brackets  may  then  be  interpreted  as  the  second  term 
in  Eq.  3.5.  A  fuller  discussion  is  given  in  Appendix  A. 

This  then  gives  us  our  set  of  rules  about  forward  versus  backward 
behavior.  As  long  as  both  indices  increase,  Eq.  3.5  may  be  used  to  ex¬ 
tend  the  sequence,  as  long  as  no  "leapfrogging"  occurs.  The  sequence  may 
be  extended  (in  either  direction)  from  the  end  point  only.  Going  backward 
in  either  index  demands  an  interpolation  formula,  now  to  be  derived. 

Let  us  start  with  and  compute  Si+n  j  and  S£  from  Eqs. 

3.8  and  3.9.  Then  compute  by  proceeding  from  either  i+n,j 

using  (3.8)  again,  or  from  i,j+m  using  (3.9).  Now  from  the  four  corner 
values  we  interpolate  for  S£+£jj+]c  (£,  <  n,  k  <  m)  at  any  interior  point 
by  the  following  algorithm: 


(S/ o) .  =  o(S/o).  .  +  8(S/o) .  .. 

vU,j+k  x,j+m 


+  ’,<s/0)i+t,,j  +  5<s/0)i^,j-h« +  eR 


(3.11) 


where  a,8,Y,S,e  are  constants  to  be  determined,  and  R  is  again  a 
completely  random  Gaussian  variate  of  unity  variance  (i.e.,  not  tied  to 
any  of  the  random  variates  used  to  compute  the  four  corner  points) . 

This  is  an  important  point — we  never  have  to  "remember"  any  previously 
used  random  variate,  in  any  step. 
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That  the  value  of  (S/o)^^ j+k  correlates  properly  with  itself 
and  the  four  corner  values  is  the  condition  that  determines  the  five 
constants.  This  leads  to  the  following  equation: 


1  =  a  +  3  +  y  +6^  +  t  +  2Pm(a8  +  y6) 


+  2Fn(ay  +  66)  +  2Fn  pm(BY  +  a6) 


.  on®  ■  ,  p^n^m  _4  k 

a  +  BP  +  yF  +  6F  P  =  F  P 


aPm  +  8  +  yfV11  +  oFn  =  fV1  k 


PF"  +  BFV  +  y  +  8p”  -  F"'1  Pk 


„n_m  .  -_n  .  „m  .  „  n-i  m-k 
aF  P  +  BF  +  yP  +  6  =  F  P 


whose  solution  is 


a  =  F£Pk(l  -  P2m_2k)(l  -  F2n  2£)/A 


B  =  FV*"k(l  -  P2k)(l  -  F2n  2£)/A 


(3.12) 


Y  =  Fn-£pk(l  -  F2£)(l  -  P2m_2k)/A 


(3.13) 


,  n- 2.  m-k / ,  _2)L  _2k. 

6  =  F  P  (1  -  F  )(1  -  P  )/A 


A  =  (1  -  F2n)(l  -  P2m) 


These  values  are  readily  seen  to  lead  to  the  correct  correlations  as 
2,k  approach  any  of  the  four  corners. 


Four-point  interpolation  may  in  some  cases  collapse  to  even  two- 
point  interpolation,  as  when  S2  1  is  computed  from  Sj  j  and  S2  2  » 
but  Eqs.  3.13  represent  the  most  general  possible  case.  The  proof  that 
this  solution  correlates  correctly  with  previously  computed  points 
(i*  <  i,  j*  <  j)  is  tedious  but  straightforward. 

We  have  thus  shown  that  two-dimensional  sequences  can  be  generated 
from  the  end  points  in  both  directions,  and  have  given  the  necessary 
equations  for  interpolating,  which  is  tantamount  to  joining  two  separate¬ 
ly  generated  sequences  as  they  near  each  other.  Note  that  from  the  form 
of  Eq.  3.5,  it  is  perfectly  appropriate  to  generate  completely  independent 
sequences  in  widely  separated  regions — the  first  term  in  (3.5)  then 
vanishes  and  only  the  random  Gaussian  is  used.  Only  when  the  sequences 
near  each  other  need  any  correlation  be  considered. 

The  previous  discussion  may  readily  be  generalized  to  more  than  two 
dimensions,  but  the  algebra  becomes  tedious  even  under  the  simplifying  as¬ 
sumptions  of  exponential  decorrelation.  Analogy  with  sequences  in  one 
and  two  dimensions  would  lead  us  to  expect  the  following  results.  Extra¬ 
polation  of  a  sequence  in  a  single  dimension  involves  only  the  end  point, 
whereas  interpolation  involves  the  two  nearest  bounding  values.  Extra¬ 
polation  in  two  dimensions  can  involve  up  to  three  previously  computed 
values  (see  Eq.  A. 4),  whereas  interpolation  in  general  demands  four  other 
points.  In  three  dimensions,  extrapolation  can  involve  up  to  seven  other 
values  (to  obtain  the  final  corner  of  a  rectangular  parallelepiped),  and 
interpolation  demands  in  general  .eight  other  points. 

We  present  below  the  case  of  extrapolation  in  three  dimensions  to 
find  the  final  value  at  the  last  corner  of  a  rectangular  parallelepiped. 
Thus  we  label  points  now  with  three  indices,  i,j,k  ,  the  last  correspond¬ 
ing  to  the  time.  Points  separated  only  in  time  are  related  by  the 
equation 


45 


(S/0) 


i,j,k+r 


+  (S/a)i,3,k+^-^  \ 


(3.14) 


where  4>  is  the  exponential  time  decorrelation.  Equation  3.5  is  replaced 
by 


(S/a)..  ,  .  =  Fnpm  <t>r(S/o).  .  , 

i+n, j+m,k+r  i,3,k 


+  ^1  -  F2n  P2m  <J,2r  R4  . 


ijk 


(3.15) 


In  Fig.  22,  let  the  axes  be  as  labeled.  Equation  3.15  is  appropriate 
only  when  no  "leapfrogging"  over  intermediate  values  occurs.  These  inter¬ 
mediate  values  (the  other  corners  shown  in  Fig.  22)  are  computed  as 
follows.  In  the  plane  k  =  k  ,  corresponding  to  a  single  instant  of  time, 
we  have  already  shown  how  to  compute  the  values  at  (i+n,j,k),  (i,j+m,k), 
and  (i+n,j+m,k)  ,  beginning  with  the  point  at  (i,j,k).  Equation  3.14  is 
then  used  to  compute  the  point  (i,j,k+r). 
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Figure  22.  Extrapolation  of  Gaussian  Sequence  in  Three  Dimensions 
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Equation  A. 4  is  then  used  to  give  (s/0)i+n  j  k+r  from  the  three 
other  corner  values  (i,j,k),  (i+n,j,k),  and  (i,j,k+r);  and  similarly  for 
the  point  at  (i,j+m,k+r)  in  terms  of  the  values  (i,j,k),  (i,j,k+r),  and 
(i,j+m,k).  Finally,  to  fill  in  the  last  corner  (i+n, j+m,k+r) ,  we  must 


write 


(S/a)..  =  A(S/°)-  •  u.  +  B(S/o) .  . 

i+n,j+m,k+r  i,j,k+r  i»3,k 


+  C(S/o)  .  +  D(S/o).  ..  ,  +  E(S/cr) .  ..  , 

i+n,3,k  i,j+m,k  i+n,j+m,k 


+  G(S/o) .  .  ..  +  H(S/o).  ..  .  ,  +  eR.  .. 

i+n,j,k+r  i,3-hn,k+r  ljk 


(3.16) 


Insuring  the  proper  correlation  with  all  other  points  leads  to  the  fol¬ 
lowing  values  for  the  coefficients  in  Eq.  3.16: 


n  m 
-F  P 


n  id  r 
F  P  $ 

m  r 
-p  $ 

— Fn$r 


(3.17) 


=  -(F2nP2tn  +  P2Vr  +  FZVr)  +  F2nP2m«p2r 
+  F2n  +  P2m  +  <}>2r 


The  formulas  for  interpolation  inside  the  rectangular  box  will  correspond¬ 
ingly  involve  the  solution  of  eighth-order  determinants. 


3.3  DIFFERENTIABILITY  OF  FIRST-ORDER  MARKOV  PROCESSES 

The  exponential  correlation  function  assumed  in  the  previous  sec- 

4 

tion  is  appropriate  to  first-order  Markov  processes,  and  is  used 
principally  for  mathematical  simplicity  in  generating  random  sequences 
in  three  dimensions.  Before  it  can  be  adopted  for  use  in  ROSCOE,  we  must 
find  a  way  of  relating  it  to  the  expected  variations  in  intensity,  while 
scanning  across  the  distribution  of  Gaussian  striations  used  in  ROSCOE, 
which  is  done  in  Sec.  3. A;  we  must  a  so  examine  more  closely  another 
problem  connected  with  differentiability. 


In  the  linear  extrapolation  of  a  sequence  of  values  of  some  quan¬ 
tity  u(x)  ,  the  simplest  is  the  first-order  sequence  for  which  u(x  +  t) 
depends  on  the  previous  value  u(x)  ,  through 


u(x  +  t) 


k  u(x)  + 


R 


(3.18) 


where  k  =  e  ■,1‘i  <1  and  R  is  a  random  Gaussian  variate.  This 
sequence,  represented  by  an  exponential  correlation  function 


3(r)  = 


2 

o 


is  pathological  in  being  continuous  but  non-differentiable.  Thus,  de- 
notfr.g  expectation  values  by  brackets  <  > 


<[u(x  +  6)  -  u(x)]2>  =  <[u(x  +  +  <[u(x)]2> 


-  2  <u(x  +  5)  u(x)> 


_  2.,  -v5. 

=  2o  (1  -  e  ) 


2 

=  2o~  vS 


(3.19) 


A8 


and,  as  6  -»■  0 


the  left-hand  side  gives 


<du/dx> 


2 

c"  being  the  variance. 


(3.20) 


The  reason  for  this  behavior  is  well  known  to  be  the  non-existence 
of  the  second  derivative  B"(0)  ,  caused  by  the  cusp  in  the  correlation 
function  3(t)  at  t  =  0  .  This  introduces  a  high-frequency  fuzz  or 
jaggedness  into  the  function  u(x)  ,  manifested  by  the  comparatively 
slow  high-frequency  fall-off  of  the  spectral  density  S(m)  : 


S(w) 


1/2- 


CO 

J *  e^*  B(t)  dr 

—  oo 


=  v/(v 


■? 

u“) 


(3.21) 


In  this  section  we  attempt  to  quantify  how  seriously  this  jagged¬ 
ness  interferes  with  the  differentiation  of  u(x)  ,  as  described  by 
Eq.  3.18,  taking  into  account  the  finite  detector  width.  In  Sec.  3.2, 
we  described  a  technique  for  generating  a  sequence  of  intensity  values 
simulating  a  striated  field  of  view  with  different  correlation  lengths 
along  and  across  the  magnetic  field,  and  also  a  finite  decorrelation 
time.  We  consider  here  only  the  one-dimensional  case,  and  let  the 
linear  field  of  view  of  the  slit  (at  the  striated  region)  be  »  ,  as¬ 
sumed  to  be  some  multiple  N  of  the  sequence  step-size  o  .  We  shall 
examine  the  character  of  the  (expectation  value  of  the)  derivative 
u'(x)  as  a  function  of  both  L  and  N  ,  where  u' (x)  is  defined  as 


u'(x)  *  (5+151  '  Uuo  +  ui  V 


+  <“-N  +  +  “0>r  >  1/2 


(3.22) 
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Letting  x  =  e  v  and  o~  =  1  ,  this  expression  may  he  shown  to 
reduce  to 


u’<*5  =  whs*  j2N  +  7i f4(N-1)x 

-  (4S  +  2)x2  +  4xN+1  +  4xN+2  -  2x2N+21 J  1/2  (3.23) 

This  is  plotted  in  Fig.  23,  and  is  seen  to  be  veiy  insensitive  to  the 
number  of  points  N  in  the  interval  &  ,  but  varies  essentially  as 

for  slit  widths  small  compared  to  the  correlation  distance  , 

in  accord  with  Eq.  3.20.  Thus  as  the  detector  field  of  view  is  narrowed 
down  to  values  small  compared  to  the  correlation  length,  the  derivative 
(as  defined  by  Eq.  3.22)  becomes  indeterminate,  as  expected. 


It  is  not  yet  obvious,  however,  that  this  spikiness  need  preclude 
the  i  .  ox  an  exponential  correlation  function  for  generating  similar 
striated  fields  of  view  even  in  this  case  (narrow  field  of  view).  It 
must  be  remembered  that  the  radiance  values  described  by  Eq.  3.18  are 
also  contaminated  by  the  other  various  sources  of  noise  (photon  noise, 
detector  noise,  etc.),  and  if  the  spectral  composition  of  the  noise  is 
much  wider  in  high  frequencies  than  that  of  the  signal,  the  derivative 
is  chiefly  due  to  the  noise.3  Thus  the  high-frequency  hash  in  Eq.  3.21 
may  itself  be  considered  as  a  kind  of  noise,  which  must  be  alleviated  in 
the  smoothing  process. 

Tne  central  reason  for  trying  to  retain  the  single-point  represen¬ 
tation  (Eq.  3.18),  with  its  admitted  bad  features  in  differentiability, 
concerns  the  algebraic  difficulties  of  extending  and  interpolating 
sequences  of  higher  order  in  two  and  three  dimensions.  In  Sec.  3.2,  we 
derived  the  equations  for  extending,  interpolating,  and  joining  different 
numerical  sequences  in  several  dimensions,  when  the  variance  of  the  com¬ 
puted  value  was  permitted  to  be  any  arbitrary  function  of  position,  the 
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statistics  were  Gaussian,  and  the  autocorrelation  function  was  exponen¬ 
tial.  Even  in  this  simplest  case  of  a  single-point  sequence,  extrapola¬ 
tion  in  three  dimensions  demanded  use  of  seven  previously  computed 
values  (to  obtain  the  final  corner  value  of  a  rectangular  parallelepiped), 
and  interpolation  between  previou: ly  computed  values  demands  the  solution 
of  eight  algebraic  equations,  in  order  to  maintain  proper  three-dimen¬ 
sional  correlation.  If  one  were  to  use  a  more  complicated  recursion 
formula  than  Eq.  3.18,  involving  more  than  the  last  previously  computed 
value,  the  difficulties  of  extrapolation  and  interpolation  would  clearly 
be  increased,  but  it  is  not  yet  clear  how  severely. 

Let  us  consider,  for  example,  an  autocorrelation  function  Gaussian 

in  shape,  since  it  is  the  cusp  at  the  origin  of  e~vlxf  that  introduces 

the  high-frequency  hash  that  leads  to  problems  in  differentiation.  Since 

a  Gaussian  correlation  function  falls  off  much  more  rapidly  than  the 

exponential,  one  might  tnink  that  in  this  case  u(x  +  t)  for  large  t 

would  be  much  less  dependent  on  the  "past"  of  the  process  than  for  the 

exponential,  which  already  corresponds  to  a  single-point  recursion  rela- 
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tion  (Eq.  3.18).  However,  the  situation  is  just  the  opposite:  for  the 

exponential  correlation,  knowledge  of  previous  values  of  u(x  -  nx) 
provides  little  information  about  u(x  +  nx)  ,  whereas  for  a  Gaussian 
autocorrelation  function,  the  value  of  the  process  in  the  arbitrary  dis¬ 
tant  future  can  be  approximated  arbitrarily  closely  by  a  linear  combina¬ 
tion  of  past  values. 


With  a  Gaussian  correlation  function 


B(t) 


2  -B2t2/4 
a  e 


(3.24) 


corresponding  to  a  spectral  density 


S(u>) 


(o' 


■n  B  e 


2  2 
-w  /e 


(3.25) 
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The  two-term  correlation  function  is  seen  to  be  regular  at  the  origin, 
and  indeed  the  quantity  u^(x)  may  be  shown  to  satisfy  the  differential 
equation 


g2  /I  2 


h  +  /2  .  , 

8 -  u2  +  U2  =  R 


(3.31) 


where  R  is  a  delta-correlated  function,  in  this  Gaussian  case  reducing 
to  white  noise.  Whether  this  process  will  actually  be  conveniently 
extendable  to  two  and  three  dimensions  has  not  been  examined. 

3.4  AUTOCORRELATION  FUNCTION  FOR  GAUSSIAN  STRIATIONS 

We  now  must  examine  how  to  approximate  the  autocorrelation  function 
appropriate  to  the  particular  (Chesnut)  distribution  of  Gaussian  stria- 
tions  used  in  ROSCOE,  by  an  exponential.  Consider  a  random  array  of 
Gaussian  rods,  with  average  axial  density  £(knf2)  ,  and  with  the  nor¬ 
malized  si.se  distribution  used  in  ROSCOE 


P(a)  = 


3  /if  0 


*  2,  2 

(a  /a)6  e~0°/a 
o 


(3.32) 


packed  in  a  slab  of  thickness  L(km)  .  Let  the  electron  density  in 
each  rod  be 


n(r)  =  n  e 
o 


2.  2 

-r  /a 


(3.33) 


_3 

so  that  the  light  emission  (cm  )  is,  assuming  for  the  moment  that  the 
rods  do  not  overlap) 


2  2 

■j  i  .  2  -2r  /a 

I(r)  =  c  n  e 
o 


(3.34) 
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Consider  the  integrated  intensity  along  a  ]ine  of  sight  (LOS)  at  P  ,  due 
to  all  the  striations  in  a  strip  dx  wide  at  a  lateral  distance  X  from 
P  .  The  probability  that  a  striation  of  size  a  is  in  dxdy  is 
P(a)  £dxdy  .  The  integrated  intensity  due  to  the  striation  of  radius  a 
is 

I(a,x)  =  c  /dy  e-<2/a2><*2  +  *2> 

-  c  a^/2  n2  e~2*2/a  (3.35) 

This  same  expression  is  true  for  all  striations,  of  size  a  in  the 
strip,  of  which  there  are  J  Ldx  .  The  autocorrelation  function  due 
to  all  striations  of  size  a  is  thus 


B(t)  =  /daP(a)  f  dxlLcWf  e-<2/a2Mx2+(x+T)2) 


/ 


2  4  9  9 

,  ,  c  no  3/2  v  3  -T  /a 

daP(a)  ~ 7 —  n  2.  L  a  e 


2  4  V  t  2 

c  n  L  to 
o  L 


°  L,v) 


(3.36) 
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Let  us  now  consider  the  case  for  which  the  detector  is  so  close  to 
the  striated  region  that  the  angular  size  of  the  striations  differs  along 

the  line  of  sight. 


Considering  again  only  striations  of  a  given  size  a  in  the  area 
dxdy  ,  the  change  in  "impact  parameter"  X  due  to  an  angular  rotation 

a  of  the  line  of  sight  is 

X'  =  X  +  a(L  +  y) 
o 


As  before  (Eq.  3.35), 
I(a,X)  = 


~(2X2/a2) 

e 


2 

cn  a 
o 


J7Fi  e 


I(a,X’) 


-(2/a2)(X2  +  2aX(Lo+y)  +  a2(LQ+y)2) 


-■■■'>  i . 


so  that 


I(a,X)  I(a.X')  =  (c2n4a2Ti)/2 

o 

e-(l/a2)(4X2  +  4aX(L+y)  +  a2(LQ+y)2) 


(3.37) 


Integrating  over  all  spatial  positions  for  the  striation  axis,  and 
over  the  size  spectrum  a  ,  gives  for  the  autocorrelation  function: 


“  >■> 

=  l  J  daP(a)  J*  dy  J  dx  I(a,X)  I(a,X')  (3.38) 


2  4  4 

c  n  ifl  ( 

O  O  V  I 

-  I  can 

3a  ( 


a(L+L)^ 


i-t 


.  -1  /aLo \  \ 

Lau 


/o  on\ 
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which  may  easily  be  seen  by  L'Hopital's  rule  to  reduce  to  (3.36)  as 
Lq  ->■  oo  ,  since  aL^  in  Eq.  3.39  corresponds  exactly  to  t  in  Eq.  3.36. 
Thus  we  have  now  formalized  the  procedure  for  integrating  along  a  line 
of  sight  that  penetrates  striated  regions  at  various  distances  from  the 
detector.  Thus,  using  a  subscript  i  to  refer  to  the  various  striated 
layers,  we  have 

B(a)  =  (it  \  a^c2)/3a 


r  4  I.  -lP<Loi+V\  .  -1  (aloi  \l 
"  |tan  \ - P - )  -  tan  j 


(3.40) 


One  of  the  objections  to  Approach  2,  the  simple  superimposition 
of  bright  structured  regions  (as  at  different  distances)  was  that  it 
apparently  locked  into  the  code  a  technique  applicable  to  detectors  at 
DSP  distances,  but  not  for  sensors  much  closer  to  the  striated  regions 
where  nearby  striations  subtend  a  much  larger  angle  than  distant  ones. 
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Equation  3.40  removes  this  objection,  and  in  Fig.  24  we  have 
plotted  Eq.  3.40  in  normalized  form  for  a  single  1000-km  thick  cloud 
at  various  distances  Lq  from  the  detector,  ranging  from  100  km  to 
essentially  infinity,  together  with  the  single  exponential.  Note  the 
similarity  of  all  the  curves;  this  agreement  is  then  a  measure  of  how 
well  the  exponential  approximation  of  Sec.  3.2  can  approximate  the  "true" 
autocorrelation  function. 

If  the  areal  density  £  of  the  rod  axes  becomes  great  enough 
that  appreciable  overlap  occurs,  correlations  between  different  rods 
must  be  considered,  since  the  intensity  of  light  emissions  is  propor¬ 
tional  to  the  square  of  the  total  electron  density,  and  thus  is  not  the 
same  as  the  sum  of  the  squares  of  the  electron  densities  appropriate 
to  the  various  rods.  In  this  case,  the  factor  1/(1  +  R^)  in  Eq.  3.36 
(R  =  t/oo)  is  replaced  by  the  expression 

1/(1  +  R 

(1  +  4F/3)2 

v  2 

where  F  =  n  /.  is  the  axial  packing  density.  For  F  -  1  ,  this 

expression  is  in  general  a  factor  of  two  or  three  higher  than  the  simple 
2 

1/(1  +  R  )  .  We  are  indebted  to  Dr.  Peter  Redmond  for  carrying  out  the 
algebra  for  this  relation. 
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Figure  24.  Autocorrelation  Function  for  Distribution  of  Gaussian  Striations 
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APPENDIX  A 


EXTRAPOLATION  AND  INTERPOLATION  OF  RANDOM  SEQUENCE 


We  have  omitted  several  steps  in  the  chain  of  proofs  necessary  to 
justify  some  of  the  statements  in  Sec.  3.2,  which  we  address  more  fully 
below. 

First,  we  must  show  the  independence  of  extrapolation  in  two 
dimensions.  Having  computed 


(S/o) 


..  .  =  Fn(S/'o).  .  +  %1  -  F2n  R. 

i+n,j  x,3  ’  i 


(A.  1) 


can  we  project  ahead  in  the  orthogonal  j  direction  from 

without  taking  account  of  S.,  .  ?  The  answer  is  in  the  affirmative, 

i+n,j 

since  if  we  try  a  solution  for  S.  of  the  form 

x,j+m 

■  ««/0)i;j  +  6(5/0^.,  +  yRj  (A.2) 


it  is  easy  to  show,  using  Eqs.  3.2  and  3.3  of  Sec.  3.2,  that 


6  =  0 


a  =  P 


(A.  3) 


so  that  the  point  at  i+n,j  does  not  influence  that  at  i,j+m  . 


Second,  we  must  modify  Eq.  3.4  of  Sec.  3.2  when  any  "leap¬ 
frogging"  occurs,  as  in  computing  the  point  at  i+n,j+m  when  we 
have  already  computed  the  other  corners  of  the  rectangle,  at  i+n,j 
and  i,j+m  ,  from  Eqs.  A.l  and  A.2.  In  this  case  we  can  easily  show 
that 
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(A.  4) 


(S/a)..  =  -Fnpm(S/a)  +  Pm(S/o) 

l+n, j+n  i,j  i+n,j 


+  Fn(S/a) .  .  +  6R. . 

1,3+m  ij 


where  R  is  a  random  Gaussian  variate  of  unity  variance,  and 


1  -  S2  .  F2n  +  P2”  -  F2np2m 


(A. 5) 


Note  the  negative  sign  of  the  first  term  on  the  right,  which  is  at 
first  sight  surprising,  since  an  exponential  decorrelation  demands 
a  positive  expectation  value  for  the  product  of  any  two  terms. 

The  answer,  of  course,  lies  in  the  positive  correlation  of  all 
the  terms  on  the  right  with  each  other,  which  must  be  considered  when 
one  takes  the  expectation  value 

<(S/o )..  ..  (S/o)  > 

i+n,j+m  1,J 

The  proof  that  (S/o)i+n, j+m  »  as  computed  in  this  manner,  cor¬ 
relates  properly  with  previously  computed  values  at  i-n",j-m”  is 
straightforward. 

The  interpolation  formula  (3.10)  has  some  interesting  features 
that  are  worthy  of  note.  If  the  interpolated  point  lies  on  an  edge, 
its  value  depends  only  on  the  values  at  the  two  corners  1  and  2 
bounding  this  edge,  and  not  at  all  on  the  opposite  corners  3  and  4. 
But  if  only  three  corners,  say  1,  3,  and  4,  have  been  computed  (and 
not  the  corner  S^)*  in  order  to  compute  the  same  point  P  cn  the 
edge  1-2,  all  three  values  at  points  1,  3,  and  4  must  be  used  in  a 
linear  combination,  as  in  Eq.  A. 4. 
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Successive  interpolation  for  points  in  a  rectangle  is  straight¬ 
forward,  but  cannot  be  done  blindly.  Having  once  computed  a  point  P 
inside  a  rectangle,  an  auxiliary  constraint  has  been  set  for  a  sub¬ 
sequent  point  P*  since  these  must  correlate  properly  with  each  other. 
Thus  Eq.  3.10  cannot  simply  be  reapplied  for  P'  .  Instead  the  two 
points  A  and  B  must  be  found  (A  from  a  linear  combination  of  2,  3,  and 
P  ,  and  B  from  a  linear  combination  of  2,  A,  and  P  ,  as  in  Eq.  A. 4), 
and  then  the  interpolation  (Eq.  3.10)  reapplied.  At  each  step  it  must 
be  shown  than  any  computed  point  correlates  properly  with  other  values 
than  those  used  in  its  determination.  See  Fig.  A.l. 

The  final  step,  of  joining  two  widely  separated  regions  computed 

wholly  independently,  is  now  seen  to  be  straightforward,  by  a  process 

of  computing  successive  corner  points  of  the  rectangle  between  the  two 
regions.  Thus  if  region  I  has  a  limiting  point  P  in  its  upper  right- 
hand  corner,  and  region  II  a  limiting  point  R  in  the  lower  left-hand 
corner,  the  other  two  corners  M  (a  linear  combination  of  P  and  R) 

and  N  (a  linear  combination  of  P,  M,  and  R)  are  then  computed  and 

intermediate  values  interpolated  as  previously  discussed.  See  Fig.  A. 2. 


A 


o> 

o 

o> 

I 


Figure  A.l.  Multiple  Interpolation  Inside  a  Rectangle 
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c 


9 


•  ACCESS  OL.  OD.  OP  1 
|  OT.  ST.  BR.  FV 
|  DATASETS  j 

|  •  INITIALIZE 
j  UPWELLING 

*  j  RADIATION 

♦ 

j  CALLUPWELL 

1 

^kflag  ol\  otherwise 

|  •  UPDATE  OBJECT  & 

\^(NEVENTS)^'>"'  - 

SENSOR  POSITIONS 

JINITIAL 

!  CALI  PLTFRM 
j  CALL  WHEREO 

! 

c 


•  CREATE  NEXT 
OPTICS  LOOK 
EVENT  AT  T  =  T  + 
FRAME  TIME 
»  PUT  AT  TOP  OF 
EVENT  LIST _ j 


•  CREATE  OPTICS 
i  PROPAGATION  EVENT 
(•  PUT  AT  TCP  OF 
|  EVENT  LiST 


n 


!  •  CREATE  SIMULATE 
1  OPTICS  EVENT 
j  #  PUT  AT  TOP  OF 
;  EVENT  LIST 


j  •  TIME  =  CLOCK  ♦  FRAME 
1  •  LOOP  ON  EVENT  LIST, 

|  IF  KTYPE  =  4,  SET 
TIME  s  TIME  E4 


•  CREATE  NEXT  LOOK 
EVENT  AT 
TIMEsTlME* 

J _ FRAME  TIME 


CREATE  LIST  OF 
n-\TH  DATASETS 
TO  TARGFTS  IN  FOV 
j  FOR  ALL  BANDS  (ILIST) 

!  •  STORE  ANGULAR  POS 
I  IN  SENSR  COORD  J 

I 


•  CREATE  OPTICAL 
!  COMPUTATION 

EVENTS  SPECIFIED  j 
!  (TRACK  HOMING. 
OlSCRIM, 


c 


T" 


)•  CREATE  REAL  LOOK 
EVENT  AT  ACO 
j  TIME.  T 
'•  KFLAGOL(NEV) 

I  =4HLOOK 


2  \ 

!  •  CREATE  SPECIAL 
j  OUTPUT  EVENT 
!  •  PUT  ON  TOP  OF 
EVENT  LIST 


NOTES:  l  FOR  ALL  NEW  PATH  DATASETS.  WE  SET  TIME  =  CLOCK 

TNEX=  -500C 

2  CURRENTLY  THERE  ARE  NO  SPECIAL  OUTPUT.  HOMING,  OR  DISCRIMINATION 
EVENTS.  THE  CODE  HAS  8EEN  STRUCTURED  FOR  ThEIR  INCLUSION  L/TER 


PROGRAM  OLOOK:  Optics  Look  Event 
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M.S6276 


CHECK  THE  STATUS  OF  EACH 
PATH  DATASET  (NSL) 

ON  HIST 


•  ACCESS  OD,  ST,  FV, 

BR,  PR 
DATASETS 

•  ACCESS  SENSOR  DATA, 
LOOK  TIME 


aFLAG  PR(NEVENTST 
\  >0? 


THIS  BLOCK  OF  THE  PROGRAM 
SETS  UP  THE  LIST  OF  PATHS 
TO  BE  INTEGRATED 


•  UPDATE 

PHENOMENOLOGY 
REGION  DATA 


CALL  INTRP 


/  SHOULD  WE  \ 
TREAT  EXTENDED 
\OBJ 


•  CREATE  LIST  OF 
PATHS  TO  OBJECTS 
FOR  ALL  BANDS  & 
PUT  ON  ILIST 
CALL  XOBJ 
CALL  XCLD 


GRID  THE 
FOV 

\  / 


•  CREATE  OR  ACCESS 
LIST  OF  PATHS  FOR 
GRID  POINTS  AT  ALL 
BANDS 

•  PUT  ON  ILIST 


0PTRP  =  3HYES 
& 

\TNEX5s  0/ 


TIME  =  TNEX 
TNEX  =  CLOCK 
KNEW=1 


PUT  PATH  ON 
(SLIST) 


'OPRP  =  3HYES'; 

\  / 


CAN  WE  \ 
INTERPOLATE:  x 
TIME<CLOCK<TNEX 

\  ^ 


•  GET  NEXT 
UPDATE  TIME 
CALL  MRCUPD 


•  CREATE  A  SECOND 
PATH  DATASET  (MSL) 
&  TRANSFER  TIME  2 
DATA  02)  TO 
TIME1  SLOT  (II) 


KNEW  =  0? 


•  CREATE  A  NEW 
TIME2  DATASET 
(12)  &  PUT  ON  (MSL) 


•  INTERPOLATE 
OPTICS  PATH 
DATA  USING  (12  & 
II  DATASETS) 

”  CALL INTRPO 


•  PUT  NEW  PATH 
DATASET  (MSL) 
ON  INTEGRATION 
LIST  (SLIST) 


PROGRAM  OPROP:  Optics  propagation  event.  Sets  up  paths  to  be  integrated 
and/or  interpolates  these  data 
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AW-56275 


•  DETERMINE  IF 
THERE  ARE  ANY 
PATHS  TO  BE 
INTEGRATED 


THIS  BLOCK  OF  THE  PROGRAM 
CLUSTERS  CLOSELY  SPACED  PATHS 


FOR  EACH 
PATH  ON  SLIST 


PATH  SL(NSL)  =  0  j  FOR  EACH  PATH 

- - '  DATASET  (NSL) 

ON  SLIST 


•  CREATE  THE  SIGHT 
PATH  CALCULATION 
EVENT  (NPT) 

AT  TIME  =  CLOCK 


FOR  EACH  PATH 
DATASET  (NSL)  ON  SLIST 


•  PUT  THIS  PATH 
DATASET  (NSL) 

ON  THE  LIST  OF 
PATHS  TO  BE 
INTEGRATED 
ORDERED  IN  TIME 
CALL  PUTORA  (ILIST 
PT(NPT).  NSL,  4) 


•  PUT  NPR  AT 
TOP  OF  EVENT 
LIST 

•  PUT  NPT  AT 
TOP  OF  EVENT 
LIST 


•  CREATE  THE  RETURN 
PROPAGATION  EVENT 
(NPR)  AT  TIME  = 
CLOCK 


TNEX= CLOCK 
OR 

KTYPE  =  6HTARGET 
OR 

PATH  *0 

NO  FOR  EACH  PATH 

i  DATASET  (MSL)  ON  SLIST 


TNEX= CLOCK 
OR 

KTYPE*4HGRI0 

OR 

PATH  *0 
XLAM  SL(NSL)* 
XLAM  SL(MSL) 
OR 

NSL  =  MSL 


NO 


PROGRAM  OPROP  (continued) 
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AN-5677J 


PROGRAM  OPROP  (continued) 
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AN-56Z73 
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NOTES 


I 


For  all  new  path  datasets,  we  set: 

TIME  =  CLOCK 
TNEX  =  -5000. 

This  part  is  computed  only  when  path  data  is  to  be  initialized. 

The  second  path  dataset  contains  the  same  data  as  the  first  for 
normal  updates.  For  initialization,  however,  a  new  021  dataset 
is  required  since  integrations  at  two  times  are  required. 

The  (II)  dataset  always  carries  the  path  data  at  the  current 
time  used  in  the  optics  simulation. 

This  part  is  computed  only  when  path  data  is  initialized.  It 
sets  up  the  (II)  data  for  interpolation. 

The  path  datasets  on  SLIST  are  ordered  such  that  initialization 
data  at  two  times  for  the  same  path  are  sequential.  A  single 
path  dataset  (NSL)  with  pointers  to  the  datasets  carrying  the 
time-varying  parameters  at  two  times  (II  and  12)  is  then  formed. 


PROGRAM  0PR0P  (concluded) 
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(  ENTER  ) 


* 

•  ACCESS  SENSOR  (OD) 

DATASET 

•  GET  SENSOR 

POSITION 

CALL  PLTFRM 

FOR  ALL  FIREBALL 

DATASETS  (NFB)  ON  FBL1S 

(  EXIT  ) 


SUBROUTINE  XOBJ: 


Generates  sight  path  datasets  to  all  fireballs,  dust 
regions,  and  beta  tubes  in  FOV 
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•  CHECK  SIZE  OF 
REGION  VERSl'S 
SENSOR  RESOLUTION 

NMAXX  =  IFIX42.  +  RM  AXi 
RESOL)- 1 

NMAXY  =  IFIX(HMAX  -  HMIN) 
RESOL)- 1 

- j 


•  CREATE  A  PATH 
OATASET  (NSL)  FOR 
EACH  SENSOR  BAND  | 

•  PUT  (NSL)  ON  PATH 
LIST  PATHLS 


SUBROUTINE  PATHS:  Routine  which  determines  number  and  position  of  paths 
to  a  specified  region 
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0119  5 -NY 


,/7\ 


CREATE  ADDITIONAL  PATH  DATASETS 
TO  REGION 


•  SET  UP  VECTOR 
NORMAL  TO  Y&VEC: 
Z  =  Y  x  VEC 


•  SET  SCALING  FACTORS: 

DIST1  =  HMAX  -  HC  -  RESOU2. 
DIST2  =  HC-HMIN  -  RESOU2. 
RC  =  (RMAX  (HC  -  60.E5)  + 
RMIN  (85. E5  -  HC)V2S.E5 


COMPUTE  POSITIONS  OF  CORNER 
POINTS  OF  REGION: 

X  s  POS  +■  DIST 1  *  VEC 
P(7.1)  =  X-RMAX‘Z 
X  =  POS-DIST2*V£C 
P(7.2)=X  +  RMIN-Z 
?-,74!=pos+rc-y 
P|7.4)  =  POS-RC"Y 
NP  =  5 


SET  UP  VECTORS  FROM  SENSOR  TO 
REGION  &  NORMAL  TO  THIS  VECTOR  & 
REGION  AXIS  (VEC). 

X  =  POS -SENSOR 
Y=XxVEC 


•  SET  UP  VECTOR  NORMAL 
TO  XSY: 

Z  =  X  x  Y 


•  SET  SCALING  FACTORS: 
ALP  =  SEPA(Z.V£C) 

AA  =  RMAX 
SB  =  HMAX - HC 
TANALP  =  TANIALP) 


DIST  =  AA*BB*  /  1  4  TANALP2 

J  aA2*  BB2  -  TANALP2 
-RESOU2. 

RMAX  =  RMAX  -  RESOLD. 


•  COMPUTE  POSITIONS  OF  FOUR  j 
POINTS  ON  REGION  PROJECTION  i 
ON  PLANE  NORMAL  TO  X:  j 

P(7.1)  =  POS4  0iST-Z 
P(7^i  =  rOS-DIST*Z 
P{7.3l  =  POS+  rmax-y 
P(7.4)=  P?=  -  RMAX-Y 
NP  =  5 


•  CREATE  ADDITIONAL  POINTS  ON 
REGION  PROJECTION  INTERMEDIATE 
TO  THOSE  ABOVE  AND: 

NP  =  7IFNMAXXW 
OR 

NMAXY>4 
NP='.3IP  NMAXXX4 
AND 

NMAXYX 


I  •  SET  UP  VECTOR 


NORMAL  TO  Y&VEC: 
Z=YxVEC 


•  SET  SCALING  FACTORS: 

DIST  =  (HMAX  —  HM1N  -  RES0L)I2. 
RMAX  =  RMAX  -  RESOU2. 


•  COMPUTE  POSITIONS  OF  CORNER 
POINTS  OF  REGION  = 

X=POS  +  DiST‘vic 

P(7.i)=x-rmax+z 

X=POS-DIST-Z  _ 

P(7.2)  =  X4RMIN  +  Z  _ 

P<7.3)=  POS  +  RMAX-Y 

P(7.4)=  POS- rmax-y 
NP  =  S 


I  •  INSURE  THAT  AT 
NO  !  LEAST  ONE  PATH 
1  LIES  IN  FOV 
IF  NOT.  DSTROY 
PATHLS 


-J 


SUBROUTINE  PATHS  (continued) 
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SUBROUTINE  MRCUPD: 


Routine,  to  determine  next  update  time  for  path 
dataset 


75 


SUBROUTINE  INTRPO:  Routine  which  is  used  to  interpolate  path  data  to 

current  time 
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AW- 56267 


j  •  ACCESS  T1  &  T? 
,  DATASETS 


nsd  >— 


WIPEOUT  OLD  L9TOT 
LIST 

Call  wiPCUT 
SET  NSAVsO 
DSPSAV  2  0 


FOR  EACH  CELL  t*> 
ON  PATH  »N$U 
I - - 


/  XLAM  SUNSL)  ;  -  I 
*  XLAMI  /  NO 
"v, “>  . ,  "  (INTEGRATE) 

GET  NEXT  (8G:  T*ES  {PUT  Oata 
DATASET)  y  AWAY) _ 

f  Q2iNI2)  =  BKTOT  BGiNBG)  j 

j  SIGMA  l^NO)  =  S!GMA  BGiNBG.  ' 

:  CORX  I2^i2i=  IK 
CORY  l*Nl2i  =  RCORE3GiN&G}  . 

:  CORTI2^i?,=  TCORE  I 

SEED  IZtNiTi  =  SEED  IliNli,  ! 


•  CREATE  THE  CELL 
LIST  FOR  THE 
NEW  PATH 
CALL  SPGEOM 


KCALC  =  (T»  > 


«CREG:K‘ 

-<  4 

'v  DlSTlOKj 


LOOP  OVER  y  -SC 

BAND  INT  tBI)  T J Af 

j  '  •"COMPUTE  PLUME 

|  :  ‘-rradiance  gu) 

I  i  CALL  £>tU«RAD 


<  •  CREATE  The  EFFECTS. i 
EMiSCAT  AND  I 

RETURN  SPCALC 
EVENTS  FOR  THIS  \ 

CELL 


/  NEXT  \  _j  •  SKIP 
\PaTh  i  INTEGRATION 


•  IF 
CJ 


"  NO 
INTEGRATE. 


IF  STRUCTURED 
OBJECT  GET 
DSPWRD  iDSPSay? 
AND  SAVE 


•  ADDI.RRAD  TO 
PATH  TOTAL 
BKGND  BliNBll »  Xj  * 
BKGND  SlMBl) 


!  •  SUM  OVER  BAND  x 
i  BKTOT  =  BKTOT  * 

|  BKGND  BHNoU 


i  •  STORE  BAND  TOTAL 
•  UP  TO  TARGET 
?  BKTOT  BG.NBG)  =  BKTOT 


PROGRAM  SPCALC: 


Sight  path  calculation  event.  Initializes  sight  path 
cell  list  and  then  creates  EFFECTS  and  EMISCAT  events 
to  process  data  cell  by  cell 


7 


. . . . . . . . . . . . 


C  ENTER 

v  .  J 


•  tNlTIAtZE  !  < 

PARAMETERS 
KPRNT  s  0 
FACT  =10.  1 

NOR  =  50 
XDELTA=  100 

_ IN  UST  a  0 _ _ 

i 

I*  COMPUTE  COSINE  OF  ; 

ZENITH:  _ 

j  „  =  SEPA(ROVEC  SPHAT) 
o  =  Tl-o  IF(.i>TU2l 
CSZEN  =  COS  («> 

''  '  I 

i  •  CREATE  THE  i 

•  BASIC  INTEGRATION  1 

LIST-  I 


SHELLS  TO  THE 
CELL  LIST- 
CALL  INSERT 


SUBROUTINE  SFGEOM:  This  subroutine  creates  the  cell  integration  list 


8 


ENTER 


D 


•  CHECK  LOOK  DIRECTION  AND 
SET  UP  VECTORS  FOR 
INTERSECTION  CALCULATIONS 
SENSVC  =  RDVEC 
OBVEC  =  2*RBODY*  SPHAT 
KFLAG  =  0  &  1=0 


SENSOR  LOOKING  UP  - 
(CONTINUE) 


•  SENSOR  LOOKING  DOWN: 

KFLAG  =  1 

or  =  ACOS(CSZEN) 

HBAR=  [RDVEC|*SINa 

\ HBAR< 

SBAR  =  VRBODY2  -  HBAR2 

|  RBODY 

,  HBAR> 

SCALE  =|RDVECI*CSZEN 

} RBODY 

SENSVC  =  RDVEC  +  SCALE*SPHAT 

OBVEC  =  RDVEC  -  SENSVC 

T 


FOR  ALL 
FBSON  FBLIS 


CHECK  FOR  NUC.  REGION  INTERSECTIONS 

I 


4HHIGH 

•  COMPUTE  EXTENT  (PL)  OF  INTERSECTION 

OF  PATH  WITH  THESE  REGIONS: 

^TYPE?^ 

CALL  EXTENT  -  FOR  FIREBALL  REGION, 

jl  3HLOW 

HA-CHEM  PANCAKE, 
CHEX-ION  LEAK. 

•  COMPUTE  EXTENT  (PL)  OF  INTERSECTION 

DEBRIS  LOSS-CONE 

OF  PATH  WITH  THESE  REGIONS: 

FOR  EACH  DEBRIS  (NDD)  $ 

CALL  EXTX  -  FOR  FIREBALL  REGION 

•  COMPUTE  EXTENT  (PL)OF  INTERSECTION 

CALL  EXTENT  -  FOR  VORTEX,  DUST 

OF  PATH  WITH  THESE  REGIONS: 

CHEILUM.  SPHERE,  WAKE 

CALL  EXTENT  ..  FOR  DEBRIS  (CONTAINED), 

CALL  SHOCKW  -  FOR  WILSON  CLOUD 

SHEATH, BETATUBE 

FOR  EACH  REGION 

♦ 

!  ^ 

n 

I _ 


SET  PARAMETERS: 

1  =  1  +  1001 
INDXWfl,  =  INDXF  FB(NFB) 
KIND(I)  = 


/10  +  KINDF-(FB) 

21 - (VORTEX) 

26 -(WAKE) 

17  -(CHM  SPHR) 

23 -(DUST)  -  , 

22 - (WILSON  CLD)'16-(CHM  PNCK) 

PUT  DISTURBED  REGION  CELL  ON  LIST  ••  CALL  INSERT 


/16-(CHX/I.L.) 

24 -(DEBRIS) 

20 -(SHEATH) 
19-(BETATUBE) 
27-(D£BRS  L.C.) 


% 


SUBROUTINE  DWELLS:  Subroutine  whic.h  creates  disturbed  region  cells  for 

specified  path 
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AM- 56264 


SUBROUTINE  DRCELLS  (continued) 
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AN-56263 


•  ACCESS  SENSOR 
DATA  SSCD.ST, 

BR.FV  DATASETS 

•  ACCESS  FOCAL  PLANE 
DATASET  LIST  JUST 


SLlST  j 


FOR  EACH  FOCAL  PLANE  ARRAY 
DATASET  (NFA)  ON  1LIST 


•  SET  UP  OET  & 

FOV  SIZES  FOR 
STORED  SENSOR 
MODELS 

•  COMPUTE  NR  &  NC 
OF  ARRAY 

~  r~ 


n 


(  EXIT 


FOR  EACH  PATH  OATASET  (Nil) 
ON  PLST  FA(NFA) 


fV  ZERO  A(lOO.lOO) 
;  &  SC(100.*,00) 

ARRAYS 


•  COMPUTE  THE  ARRAY 
INDICES  OF  POINT 
AND  SET  ARRAY 
VALUE  TO  PATH  ? 

VALUE _ ___ 

ACKK.LL)  =  Hlt(NII) 


•  PERFORM  SENSOR 

PROCESSING _ 

•  SET  UP  SPIRE  & 

RETURN  OPSIM  EVENTS 
AT  TIME  =  CLOCK 


/next\ 

VS/ 


FOR  EACH  PATH  OATASET  (Nit) 
ON  PLST  FA(NFA) 


GENERATE  LOCAL 
REGION  ABOUT  TARGET 
IN  FOCAL  PLANE  ARRAY 


PERFORM  SENSOR 

PROCESSING _ 

SET  UP  SPIRE  & 

RETURN  OPSIM 
EVENTS  AT  TIME  =  CLOCK 


GENERATE  REPRESENTATION  j 
OF  THIS  OBJECT  IN 
!  ARRAY  USING  PTSLS 
|  CALL  ARRAY 


f next\ 
\  NFA  j 


•  PERFORM  SENSOR 
PROCESSING 

•  SET  UP  SPIRE  4  RETURN 
OPSIM  EVENTS 


PROGRAM  OPSIM: 


Optics  simulation  event.  Generates  focal  plane  array  and 
calls  processing  routines 
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PROGRAM  OPSIM  (continued) 


82 


FOR  EACH  PATH  DATASET  (Nil) 
ON  H  IST _ 


FOR  EACH  CELL  (I.J)  IN  ARRAY 


(  Ex7t~) 


•THESE  COORDINATES  ARE  IN  THE 
ROTATED  FRAME  (ALIGNED  WITH 
THE  FIELD)  FPCOft 


SUBROUTINE  ARRAY :  Constructs  the  focal  plane  intensity  array,  given  path 
data  describing  target  prints,  objects,  and  grid 


S3 


•  SET  UP  INDICES  IN  FPCOR  COORDINATES 
FOR  REGION- 

NX  =  IF|X(AM  AX  -  AMIN  V{2  X  01  AM  +  1  )  x  2  +  1 
NY  =  IFIXIEMAX  -  EMINV12  X  DIAM  +  1  )  X  2  *  1 


FOR  EACH  CELL  (KK.  LL) 
IN  FPCOR  COORDINATES 


•  CONVERT  TO  SENSOR 
COORDINATES- 

XXX(1)  =  RA(1)  +  FLOAT(KK  -  (NX  +  1V2)‘ DIAM 
XXX(2)  =  RE(1)  +  FLOAT(LL  -  (NY  +  I)f2)-0IAM  | 
CALL  FPCOR  (XXX.  RAE.  .. ) 


NEXT  \  NO  . 
CELL 


CELL  IN  FOV. 
RAE(2)>A2M/2 
OR 

RAE(3)>ELV/2 


DOES  PATH  THROUGH 
THIS  CELL  INTERSECT 
REGION  (PL* OP 
CALL  EXTENT 
OR 

CALL  CLOUD  2 


LOOP  OVER  PATH 
DATASETS  (Nil)  ON 
OBJECT  LIST  ILIST 
TO  FIND  THREE  PATHS 
CLOSEST  TO  THIS  CELL 
XJ(K).  RA(K).  RE(K) 

K  =  1,3 


CALCULATE  DIFFERENCES 

XJ1=XJ(1)  SlGl=SIGMA(l)  CORX1  s 

XJ2  =  XJ(2)  SIG2  =  SIGMA(2) 

XJ3  =  XJ(3)  SIG3  =  SIGM  A(3) 

X24  =  RA(2|  -  XXX(l)  Y41  =  XXX(2)  -  RE(I) 


STRUCT  ^ 
=  3HYES"' 


I  •  INTERPOLATE  FOR  INTENSITY  VARIATION.  AND  i 

SPATIAL  CORRELATIONS 

ASIGMA  =  (SIG1  X  FAC1  +  SIG2-FAC2  +  SIG3'FAC3)/FAC4  | 
ACORX  =  ICORXI-FACt  +  CORX2-FAC2 ♦  CORX3*FAC3L'FAC4  l 
ACORY  =(CORY1’FACl  +  CORY2*FAC2+ CORY3*FAC3)/FAC4  I 
F  =  EXP(  -  DIAM/ACORX)  I 

P=EXP(-  DIAM, 'ACORY) 


•  COMPUTE  RANDOM  SAMPLES 
F  =0  IF  SC(KK  -  1.LL)  =  0 
P  =  0  IF  SC(KK.LL-1)  =  0 
S001  =  SC1KK  -  1.LL  -  1) 

S101=SC(KK.  LL-1) 

SOU  =SC(KK-  l.LL) 

EPS  =  0 

EPS1  =  (1  -  (F2  +  P2  +  T2  *  F2P2T2  -  (F2P2  +  F2T2  *  P2T2))) 
EPS=  \EP$i  IF  EPS1M) 

EPS2  =  (1  -  F2  -  P2  +  F2P2) 

EPS  =  EPS  +  T ' N  EPS2  JF  EPS2>0. 

SC(KK,  LL)  =  -  F'P'SOOI  +  P'LSIOl  *  F'SOII 
+  EPS*RNV(0 ,1 ) 


•  COMPUTE  ARRAY  INTENSITY  AT  THIS  CELL 
AOLD  =  A(I.J) 

A(I.J)  =  AMEAN  *  ASIGMA’SCIKK,  LL) 

A(I.J)  =  AMAX1(AOLD.  A(I.J)) 


Y23  =  RE(2)  -  RE(3)  X2l  =  RA(2)  -  RA(I) 

X32  =  RA(3)  -  RA(2)  Y2«  =  RE(2)  -  RE(1) 

Y24  =  RE(2)  -  XXX(2)  FAC1  =  X24*Y23  +  X32'  Y24 

X41  =  XXX(I)  -  RA(1)  FAC2  =  X41  *Y!3  +  X3P  Y41 
Y13  =  RE(1)  -  RE(3)  FAC3  =  X41 -Y21  -  X2I -Y41 

X31  =  RA(3) -  RAC!,.  FAC4  =  X21 -Y23  +  X32’Y2; 


Y13=RE(1)-RE(3) 
X31  =  RA(3)-RA(’, 


SUBROUTINE  ARRAY  (continued) 
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•  ACCESS  SENSOR  DATASETS: 
OD,  ST,  FV,  BR 

•  ACCESS  FOCAL  PLANE  ARRAY 
AA(!,J) 

.  ;  i  - 

•  CONVERT  NORMALIZED  ARRAY 
INTENSITIES  TO  RADIANCE 

IN  (WATTS/CM2/SR) 

|  A(I,J)  =  A(I,J)*NEFD/DIAM2 
FOR  ALL  1,J 

1  j 

•  ACCOUNT  FOR  BLURRING 
OF  ARRAY  DURING  SCAN: 
CALL  ABLUR 


•  GET  INDEX  OF 
NEXT  TARGET 
ON  TARGET 
LIST 

(DSP  PC(NEVENTS) 


•  CREATE  THE  OUTPUT 
DATASET  (NOM) 

•  ZERO  THE  OUTPUT  ARRAY 
A(l,J)  =  0. 


•  FOR  EACH  PROCESSING  BLK 
ON  SPIRE  LIST  CALL  PROPER 
ROUTINE  TO  PROCESS  IT: 
CALL  SPSUB1 
OR  CALL  SPSUB2 


TYPE  ^ 
:  TARGET?^. 


(  NEXT  1 
i  TARGET i 


LAST 

TARGET 


CREATE  THE  TARGET  MEAS.: 
AZMES  OM(NOM)  =  SUMXO 
ELMES  OM(NOM)  =  SUMYO 
RAD  OM(NCM)  =  ARAD"  NEFD 
ETC. 


PROGRAM  SPIRE:  This  program  calls  the  subroutines  to  scan  the  focal  plane 
and  process  the  data.  Also  produces  ome  output  datasets 
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MODEL 

TYPE 


OTHER 


HERE  WE  HAVE  LEFT  SPACE  FOR 
STORED  SCAN/PROCLSSING  BLOCK 
MODELS 


FOR  EACH  BLOCK  (NOl) 
ON  SPIRE  LIST  LHV 


r  • 

•  SET:  BLKS  =  BLKS  +  1 

•  READ  DATA  FROM  j 

SET: 

1BLK  =  IBLK  +  1 

U'ARRAY  INTO  BLKS  j  _ _ __ 

NDET  r  NDET  ST(NST) 

•  GET  NUMBER  OF  WORDS  IN 

ARHat  j  /NEXt\ 

BLURO  =  SQRT(BLK(2.1)**2 

THIS  BLOCK: 

BLK(I.IBLK)=  !-n  NOl  J 

+  BLK(3.1)*‘2R 

IDX  =  IFIX  (Q(NOl)) 

OQ(NOl.l)  j  V - / 

SCS  =  BLURO 

NMBR  =  ND0(1DX) 

FORIsl.NMBR 

f 

OUTFLG  =  2. 

SUBROUTINE  INSENS:  Initializes  SPIRE  scan/processing  computation  blocks 


•  COMPUTE  BRANCH 
NO.  &  INDICES: 

BRI  =  INT  (BRANCH +  .1) 

LI  =  1 

L2  =  (SCAN  LENGTH/SAMP) 


COMPUTE  SCAN  RATE  (SCRA) 
SAMPLE  TIME  (SAMPT),  NORM. 
FACTOR  (DNORM): 
cpoa-  SCAN  LENGTH 
SCAN  PERIOD 
SAMPT  =  SAMP/SCRA 
DN0RM  =  6  28*BLUR02* 
SQRT(DNEB-2*SAMPT)/ 

DNEF  I 


•SCAN  EQUATIONS 
(SCAN1:  CIRCULAR  SCAN) 

X  =  ARRX/SIN(TH)  x  SINfT/PH  x  SL  +  TH) 

Y  =  YCOORD  -  SIGN(ABS(YCOORD  -  ARRYyCOSfTH) 
x  COS(T/PH  x  SL  +  TH).  YCOORD) 

(SCAN2:  LINEAR  SCAN) 

X  =  ARRX  +  SCRA  x  T 

Y  =  ARRY 

(SCAN3:  DSP  CIRCULAR  SCAN) 

X  =  ABS(YCQORD  -  ARRY)  x  SIN)  62832  x  T) 

Y  =  YCOORD  -  SIGN(ABS(YCOORD  -  ARRY)  x  COS 
( 62832 *T),  YCOORD) 

T  =  .62832 "T 
(SCAN4-  DUMMY) 

X  =  Y  =  T  =  0. 


“L“  LOOP  OVER  NO. 
OF  DETECTORS  (NDET) 


SET  BRANCH  NO 
&  INDICES: 


•J"  LOOP  (FROM  LI  TO  L2) 
FOR  SCAN  INTERVAL 


•  COMPUTE  SCAN  POSITION  FOR  j 
THIS  SAMPLE: 

X  =  ARRX(J)  &  Y  =  ARRY(J) 

Ts(J-  100)*SAMPT 

•  ACCCESS  PROPER  SCAN  FUNCTION 
DEPENDING  ON  SCNFLG: 

SCAN1  FOR  SCNFLG  =  1.  * 

SCAN2  FOR  SCNFLG  =  2. 

SCAN3  FOR  SCNFLG  =  3. 

SCAN4  FOR  SCNFLG  =  4 


•  CONVERT  COORD  TO  ARRAY. 

X1S  =  X1  +  XFV02 

Y1S  =  Y1  +  YFV02 

•  GET  ARRAY  INDICES  AT  THIS 
POSITION: 

NX  =  INT(X1S/SCS  +  0.5)+  1 
NY  =  INT(Y1SlSCS  +  0.5)  +  1 

•  COMPUTE  NORMALIZED  SIGNAL; 
A(J.BR)  =  SC(NX,  N Y)  +  DNORM 


I  NEXTL. 

\  i  1 


BR  =  BRI  +  (L  -  1)' 
BRANCH  INTRVL 
LA(BR)  =  LI 
LB(BR)  =  L2 


•  CREATE  A  SAMPLE  WAVEFORM  DATASET  (NSW) 
AND  PUT  ON  OUTPUT  (SAM PL  OM(NOM))  LIST: 
CALL  CREATE(NMBR  SW.  NSW) 

TIME  SW(NSW)  =  TOUT  OM(NOM)+  FLOAT(J  -  100)* 

SAMPT 

DET  SW(NSW)=  FLOAT(BR) 

KSW(NSW)  =  J 

BAND  SW(NSW)=  BAND  OM(NOM) 

A Z  SW(NSW)  =  XI 
ELSW(NSW)  =  Y1 
RAD  SW(NSW)  =  SC(Na.  NY) 

RADN  SW(NSW)  =  A(J,BR) 

*  CALL  PUTBOT  (SAMPL  OM(NOM).NSW) 


SUBROUTINE  SPSUB1  — 

(4)  Scan  Block — assuming  square  detector 
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SUBROUTINE  SPSUB1  — 
(5)  Detector  Noise  Block 
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AN- 56 2 56 


ENTER 


•  SET  OUTPUT  VARIABLES: 
J  =  K  SW(NSW) 

RADS  SW(NSW)  =  A(J,BR) 


SUBROUTINE  SPSUB2  — 

(6)  Gaussian  Smoothing  Block 
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AN-56252 


•  ACCESS  DAT AS£t3  j 
NOB.  NOD.  NOP.  NBR.  j 


/'"ktrak  TS 
(NEVENTS) 


PROGRAM  OTRACK  -  THIS  PROGRAM  SIMULATES  THE  TRACKING  OF 
AN  OBJECT  BY  AN  OPTICAL  SENSOR.  THESE 
OPTIONS  ARE  ALLOWED: 

©  AUTONOMOUS  SENSOR  TRACK  -  SENSOR 
IS  GIVEN  SOME  RANGE  EST.  FROM 
ANOTHER  SOURCE  &  TRACK  IS  INITIATED 
AUTONOMOUSLY  (IACQ  BR(NBR)  =  3HYES  & 
j  SNET  BO(IBASE)  =  2HNO) 

I  ©  NETTED  SENSORS  -  SENSO.' MAY  PERFORM 

*■'  ACQ.  &  MEAS.  ARE  NETTED 

(IACQ  =  3HYES.  SNET  =  3HYES) 

@  NETTED  SENSOR  NO  ACQ  -OPTICAL SENSOR 
CANNOT  ACQ  (IACQ  =  2HNQ.  SHET  =  3HYES) 


IACQ  BR(NBR)  > 


AUTONOMOUS  SENSOR  TRACK  ••  1ST  LOOK 

•  GET  MEAS  DATASET  (NOM)  FOR 
OBJECT  CORRESP  TO  THIS  EVENT 

•  PRODUCE  RANGE  MEAS  STAKE 
ELV  S  AZM  MEAS  FROM  (NOM) 

•  CONVERT  TO  XYZ  COORD 
-  CALL  OTRAN 

•  CREATE  A  TEMPORY  TRACK  FILE  OF 
A  WORDS  TO  HOLD  TIME.  POSIT 
WITH  DSP  WORD  =  KOSPW 

•  FIND  NEXT  LOOK  EVENT  ON  EVENT 
LIST  S  SET 

KTRAK  OL(NVNT)  =  KDSPW 


SUBSEQUENT  LOOKS 

•  COMPUTE  LENGTH  OF  TRACK  FILE  S  SET 
KOUNT  FLAG  = 

i  1  N  WDS  -  A  ..  ONE  LOOK 

(2  NWDS-  3 -TWO  LOOKS 
13  NWDS  =12-  THREE  LOOKS 
U  NWDS>12  -  TRACK  FilF  ESTABLISHED 

•  GET  MEAS.  DATASET  (NOM)  FOR  OBJECT 
CORRESP.  TO  THIS  EVENT 

•  PRODUCE  RANGE  MEAS.  FOR  KOUNTC3 
OTHERWISE  USE  EST.  RANGE  S  GET 
ELV  S  A2M  MEAS. 

•  CONVERT  TO  XYZ  COORD 


2ND  S  3RD  LOOKS 

•  EXPAND  TRACK  FILE  TO  INCLUDE 
THIS  MEAS.  -  NWDS  =  NWDS  +  A 

•  FIND  NEXT  LOOK  EVENT  ON  EVENT 
LIST  5  SET  - 

KTRAK  OL(NVNT)  =  KTRAK  TS(N EVENTS) 


;  •  SET  UP  SlGMAT  MATRIX  A  SIGMA  ARRAYS: 
j  SIGMAT(1.1)  =  RERRTS<NEVENTS)**2 

■  SIGMA(1)  =  SQRT(S:GMAT(1,D) 

I  SIGMAT(J_J)  =  FIXESE(NSE  +  J-2)"2 
|  TSNESE(NSE+J~2)"2SNR 

■  SIGMA(J) = SQRT(SIGMAT(JJ)) 


PROGRAM  OTRACK:  Simulates  the  tracking  of  an  object  by  an  optical  sensor 
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PROGRAM  OTRACK  (continued) 
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